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Applications  to  Reliability.  (August  1976) 

Jon  0.  Epperson,  B.S.,  University  of  New  Mexico; 

M.A. , University  of  Missouri 
Chairman  of  Advisory  Coaxlttee:  Dr.  J.H.  Metis 

The  cuaulant  generating  functions  of  some  multi-compartment 
stochastic  models  having  time-dependent  transition  rates  are  derived. 
Using  these,  the  first  and  second  moments  of  the  stochastic  distribu- 
tions are  found.  Under  certain  assumptions  concerning  the  Initial 
conditions  the  stochastic  distributions  are  specifically  Identified. 

The  models  Investigated  are  the  n-compartment  mixed  catenary-mammillary 
nodal  and  the  n-conpartment  general  Irreversible  model.  These  deriva- 
tions are  then  used  to  analyze  some  reliability  systems.  In  particular, 
the  catenary  portion  of  the  mixed  model  Is  used  to  obtain  stochastic 
analyses  of  a standby  redundant  system.  The  mixed  model  la  used  for  a 
system  wherein  a component  Is  subject  to  competlng-rlsks  with  standby 
co^Kments.  The  general  model  Is  employed  to  evaluate  a mission 
reliability  with  a hierarchical  failure  structure. 
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1.  INTRODUCTION 
1.1  Preliminaries 

Compartmental  analysis  is  a relatively  new  branch  of  mathematical 
modeling  which  has  undergone  considerable  growth  in  recent  years. 
Contributions  to  this  theory  have  come  from  diverse  scientific 
disciplines  which  attests  to  the  generality  and  wide  applicability  of 
the  method.  In  the  general  compartmental  modeling  problem,  one  wishes 
to  describe  'che  movement  of  a substance  through  a system.  This  may, 
for  example,  be  the  movement  of  a chemical  through  an  organism,  a 
person  through  an  organization  or  perhaps  an  organism  through  an  eco- 
system. Whatever  the  situation,  it  is  assumed  the  system  under  con- 
sideration can  be  divided  into  homogeneous  compartments.  The  rate  at 
which  the  substance  leaves  a conq>artment  is  usually  assumed  to  be 
proportional  to  the  amount  of  substance  in  that  compartment.  Other 
mathematical  formulations  for  compartmental  transitions  have  been  in- 
vestigated and  are  suomarized  in  Section  1.3.  Thinking  in  terms  of  a 
unit  or  particle  of  this  substance,  it  is  further  assumed  that  one 
unit  in  a compartment  acts  Independently  of  every  other  unit  in  that 
compartment. 

It  is  clear  these  assumptions  greatly  restrict  the  problems 
which  can  be  modeled  in  this  way.  However  the  class  of  problems 
which  fits  these  assumptions  is  still  quite  broad.  Sheppard  [1962], 
Rescigno  and  Segre  [1965],  Atkins  [1969]  and  Jacquez  [1972]  discuss  a 


Citations  will  follow  the  format  of  Biometrics . 
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wide  variety  of  such  problems. 

In  classical  compartmental  analysis  the  transitions  from  one  com- 
partment to  another  are  assumed  to  be  deterministic  In  nature.  This 
assumption  leads  to  a system  of  linear,  first-order  differential 
equations  whose  solution  yields  the  amount  of  substance  in  a given 
compartment  at  any  time.  Many  systems  have  been  adequately  modeled  In 
this  way.  However,  In  some  cases  more  realism  can  be  put  Into  the 
model  by  assuming  a stochastic  or  probabilistic  behavior  for  the  com- 
partmental transitions.  Matls  [1970]  states  this  advantage  of  the 
probabilistic  approach  In  the  following  way; 

...a  compartmental  model  can  only  claim  to  represent 
an  'approximate'  theoretical  background  for  the 
biological  observed  phenomena.  In  that  It  employs  an 
abstraction  of  'transfer'  of  particles  without 
specifying  a detailed  causative  theory  responsible 
for  the  transfer  mechanisms.  It  Is  also  because  of 
this  approximate  feature  that  a stochastic  compartmental 
theory  Is  more  realistic  In  that  a detailed  causative 
mechuilsm  %«hlch  Is  lacking  In  the  deterministic  model 
Is  replaced  by  a stochastic  model. 

The  stochastic  model  not  only  has  this  added  realism  but  In  some 
cases  one  can  recover  the  deterministic  model  In  that  the  expectation 
of  the  stochastic  model  Is  the  same  as  the  deterministic  solution. 

Such  a circumstance  Is  discussed  In  Matls  [1970].  This  is  not  to  say 
stochastic  modeling  will  or  should  replace  deterministic  modeling. 

Matls  [1976]  argues  that  while  some  problems  are  adequately  modeled 
deterministically,  other  problems  currently  approached  deterministically 
could  be  iBodeled  better  stochastically. 

Uppulurl  and  Bernard  [1967]  and  Grandljan  and  Bargnar  [1972] 
discuss  stochastic  models  whose  expectations  are  not  the  same  as  their 
corresponding  deterministic  models. 
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1.2  ApplicatloDB 


la  this  section  eusqples  will  be  presented  to  Illustrate  the 
variety  of  probleas  to  which  the  coapartaental  technique  has  been 
applied.  These  exaaples  are  divided  Into  deteralnlatlc  aodels  and 
stochastic  models.  It  la  convenient  to  Illustrate  a compartmental 
model  using  a schematic  wherein  boxes  represent  compartments  and  arrows 
represent  Intercompartmental  transitions.  The  transition  rate  from 
compartment  j to  compartment  1 Is  denoted 


1.2.1  Applications  using  Deterministic  Models 

Atkins  [1969]  gives  a comprehensive  account  of  the  use  of  determin- 
istic methods  In  biological  systems.  One  of  his  examples  (p.  69}  Is  the 
two-compartment  modal  depicted  in  Figure  1.  This  was  used  to  model  the 
movement  of  water  within  a rabbit.  Here  the  primary  question  Is  not 
how  much  wster  Is  In  each  compartment  but  rather  the  kinetics  of  the 


Figure  1 

Two-Compartment  Model  for  Water  Movement  Within  a Rabbit 
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water  within  the  system.  To  do  this  the  water  must  be  labeled  and  the 
movement  of  the  labeled  std>stance,  called  tracer.  Is  observed. 

In  this  study  inulln  was  used  as  a tracer  substance.  A single 
Intravenous  Injection  of  Inulln  was  given  and  the  concentration  of 
inulln  within  the  plasma  determined  at  subsequent  time  periods.  These 
data  then  allowed  for  an  estimation  of  the  exchange  rates  ^21^ 

of  water  between  the  plasma  and  extracellular  fluid  and  the  rate  of 
excretion  of  water  as  urine  . 

This  example  is  typical  of  a situation  which  often  occurs.  That 
Is,  one  compartment,  the  plasma.  Is  readily  accessible  to  the  researcher 
whereas  the  other  compartaient,  the  extracelltilar  fluid,  is  not  so  acces- 
sible. The  compartmental  approach  often  allows  estimates  to  be  made  on 
regions  of  a system  where  observations  cannot  be  made  directly. 

The  compartmental  approach  was  used  by  Herbst  [1963]  to  model 
organizational  connltment.  He  assumed  the  decision  structure  shown  In 
Figure  2.  Using  time  series  data  on  the  departures  of  Individuals  from 
firms,  he  was  able  to  estimate  the  transition  rates  as  well  as  the 

number  of  persons  In  each  compartment  at  time  t.  Note  that  this  example 
also  Includes  compartments  not  readily  observable  but  irhose  content  Is 
Invaluable  to  a firm's  personnel  recruitment  planning.  This  Is  an 
example  In  which  the  deterministic  approach  provided  an  excellent  fit  of 
the  data. 


5 


Figure  2 

Organizational  Coanitmenc  Model  as  a Conpartaental  System 
1.2.2  Applications  using  Stochastic  Models 

Using  indigestible,  plastic  beads  as  tracers.  Metis  [1970]  was 
able  to  estimate  the  parameters  for  the  ruminant  gastrointestinal  tract 
illustrated  In  Figure  3.  (In  addition  he  evaluated  the  fit  of  this 
model  and  found  a two-coaq>artment  model  to  be  superior.) 

At  time  t ■ 0,  4,000  of  the  beads  were  placed  into  a sheep's  rumen. 
The  feces  were  collected  at  fixed  time  intervals  and  bead  counts  made. 
Using  these  time  series  data  the  transition  rates  were  estimated  by  a 
method  which  made  use  of  the  stochastic  conpartmental  distribution 
and  the  serial  correlation  of  the  data. 
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Feces 


Figure  3 

BuBioant  Gastrolntesclnel  Tract  as  a Cotq>artaental  Model 

For  a final  exaaple  of  the  use  of  stochastic  conpartaental  model- 
ing, consider  the  compartmental  arrangement  of  Figure  4.  This  Is  known 
as  a niiwminsry  model  and  has  an  abundance  of  applications.  Matls 
at  al.  [1974]  consider  the  central  compartment  to  be  the  blood  stream 
and  tha  peripheral  compartsttnts  to  be  tissue  groups,  k particle  In 
the  blood  stream  can  pass  to  ons  of  the  k tissue  groups  with  rate  a 
or  leave  the  system  completely  with  rate  0.  Once  a particle  has  entered 
a tissue  group  it  can  leave  the  system  with  rata  y. 
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Figure  4 

A Kaomlllary  Conpartaental  Model  with  k Peripheral  Compartments 


Using  stochastic  theory  they  determine  lower  bounds  on  the 
probability  of  having  a given  number  of  particles  present  In  a given 
tissue  group.  This  application  has  Implications  in  cancer  research 
under  the  assumption  that  a carcinogen  present  In  the  blood  stream 
Induces  cancer  In  a tissue  group  If  a critical  threshold  level  of  the 
substance  is  reached  In  that  tissue  group.  This  solution  requires 
ths  use  of  the  stochastic  theory  and  cannot  be  solved  with  the 
deterministic  model  alone. 
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1 . 3 Literature  Review 


This  dissertation  Is  prloarlly  concerned  with  the  derivation  of 
some  particular  stochastic  models  and  their  applications.  The  deter- 
ministic coimterparts  to  these  models  will  play  a role  In  the  appli- 
cations but  are  derived  by  distribution  theoretic  considerations  rather 
than  methods  found  In  the  literature  of  deterministic  compartmental 
analysis.  This  literature  is  of  Interest  here  due  to  Its  place  In  the 
evolution  of  the  compartmental  theory. 

Atkins  [1969]  was  previously  referenced  for  his  account  of  deter- 
ministic methods  in  biological  systems.  Other  references  for  the 
deterministic  method  are  Sheppard  [1962],  Resclgno  and  Segre  [1965] 
and  Jacquez  [1972]. 

Bartholomay  [1958]  was  one  of  the  first  to  suggest  the  stochastic 
approach  and  solved  the  one-compartment  model.  Since  the  Initial  work 
of  Bartholomay,  the  theoretical  evolution  of  stochastic  compartmental 
analysis  has  proceeded  In  two  general  categories.  On  the  one  hand 
generalizations  have  been  made  to  the  mathematical  formulation  of  the 
transition  rates  and  on  the  other  hand  modifications  and  extensions 
have  been  made  to  the  compartmental  structure  Itself. 

Initially,  transition  rates  were  taken  to  be  constants.  These 
have  since  been  modified  In  a number  of  Independent  ways.  Uppulurl  and 
Bernard  [1967]  Introduce  the  Idea  of  a random  component  between  compart- 
ments and  Chlang  [1968]  Incorporates  time-dependent  Input  rates.  A 
two-compartment  system  having  all  transition  rates  time-dependent  is 
solved  by  Purl  [1968].  However,  his  solution  Involves  a series 
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representation  for  the  probability  generating  function.  Concentration- 
dependent  transition  rates  are  Introduced  by  Bellman  [1970].  Sykes 
[1969]  and  Soong  [1971,  1972]  treat  the  transition  rates  as  random 
variables . 

The  one-compartment  model  Is  generalized  by  Thakur  [1972]  such 
that  the  compartmental  Input  Is  time-dependent  and  In  addition,  the 
usual  assumption  that  the  Initial  number  of  particles  be  known  Is 
dropped.  Instead  the  more  general  asstmptlon  Is  made  that  their 
distribution  Is  known.  One  generalization  which  has  received  less 
attention  than  those  applied  to  the  transition  rates  Is  the  Introduc- 
tion of  age-dependency  to  compartments.  Blaster's  [1956]  model  for 
passage  of  particles  through  the  gastrointestinal  tract  of  ruminants 
Is  generalized  by  Hatls  [1972]  to  Include  gamma  age-dependency  In  the 
first  compartment.  Purdue  [1974]  Investigates  one  and  two-compartment 
models  In  which  particles  present  at  time  zero  have  different  lifetime 
distributions  from  those  arriving  later.  A.  queuing  theoretic  approach 
Is  used  In  this  analysis. 

Results  which  are  advancements  in  compartmental  structure  of 
course  begin  with  the  solutions  to  the  one  and  two-compartment  systems 
which  assume  constant  transition  rates.  Katls  and  Carter  [1969]  find 
the  first  and  second  moments  of  the  two-compartment  system.  Katls  and 

I] 

Hartley  [1971]  solve  the  general  n-compartment  model  with  constant 

transition  rates  and  propose  a method  for  estimating  the  parameters  t 

I based  on  Hartley's  [1961]  modified  Causs-Newton  method  for  fitting  non- 

I linear  regression  functions  by  least  squares.  In  this  model  a particle 

' In  one  compartment  can  transfer  to  any  other  compartment.  Hence  the 

I 
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flow  is  completely  reversible.  Cardenas  and  Matis  [1974]  solve  the 
Irreversible  n-compartment  catenary  system  and  the  irreversible  n-com- 
partment  mammillary  system  assuming  time-dependent  input,  output  and 
transition  rates.  These  results  are  in  turn  used  in  Cardenas  and  Katls 
[1975]  to  solve  the  general  reversible  two-compartment  model  with  time- 
dependent  input,  output  and  transition  rates.  A special  class  of  n- 
compartment  systems  is  investigated  by  Cardenas  and  Hatls  [1975]  wherein 
the  structure  is  general  and  reversible  but  the  transition  rates  are 
some  fliultlple  of  a time-dependent  function. 

Matis  [1976]  presents  an  overview  of  deterministic  and  stochastic 
techniques.  In  doing  so  he  discusses  three  general  steps  in  applying 
the  coaq>artmental  method  to  specific  systems.  The  first  step  is  the 
development  of  a plausible  and  relevant  compartmental  model  for  the 
system  of  interest.  The  second  step  is  the  derivation  of  the  mathe- 
matical solution  to  the  proposed  model  and  the  final  step  consists  of 
obtaining  the  optimal  statistical  design  for  data  acquisition  and  then 
estimating  the  parameters  of  the  model  from  these  data. 

The  estimation  phase  is  by  far  the  most  difficult  and  the  least 
developed  of  these  three  steps.  See  Kodell  [1974]  for  a comprehensive 
review  of  estimation  procedures  in  stochastic  modeling. 

1.4  Assumptions  and  Notation 

Let  X^(t)  denote  a stochastic  variable  indicating  the  number  of 
units  in  the  i^^  compartment  of  a system  at  time  t.  The  system  will 
consist  of  a fixed  number  of  compartments,  n,  and  compartment  zero  will 
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Indicate  the  system  exterior.  The  transition  rate  from  the  compart- 
ment to  the  i^*^  compartment  at  time  t will  be  denoted  X^j(t).  The 
change  in  the  number  of  units  in  compartment  i during  the  time  interval 
(t,  t+At)  is  AX^(t),  i.e. , 


AX^(t)  - X^(t+At)  - X^(t) 

for  i - 1,2,  . . . , n. 


It  is  assumed  units  transfer  in  any  of  three  ways;  a unit  may  enter 
a compartment  from  the  system  exterior,  move  from  one  compartment  to 
another  compartment  or  exit  a compartment  to  the  exterior  of  the  system. 
These  transfers  have  various  names  depending  on  the  discipline.  For 
example,  the  ecologist  would  speak  of  these  as  immigration,  migration 
and  death,  or  the  biologist  as  ingestion,  circulation  and  excretion. 

The  probability  of  these  three  events  occurring  will  now  be  defined  as 
follows : 


Prob  {a  single  unit  enters  compartment  J 

from  outside  the  system  In  the  interval 
(t,  t+At)}  - 

XjQ(t)At+o(At)  j - 1,2,  ...,  n, 

Prob  {a  single  unit  in  compartment  J moves  to 

compartment  1 in  the  Interval  (t,  t+At) } - 
Xj (t) (t) At  + o(At)  for  all  i and  J, 
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and 

Prob  {a  single  unit  In  compartment  j leaves  the 
system  In  the  Interval  (t,  t+At) } •• 

Xj(t)XQj(t)At  + o(At)  for  all  j. 

Note  that  the  probability  of  two  or  more  events  occurring  during  the 
Interval  (t,  t+At)  Is  o(At) . 

The  following  assumptions  are  made  throughout  the  dissertation: 

1)  The  units  In  a compartmental  system  act 
Independent  of  each  other. 

11)  Units  entering  the  system  from  the  outside  do 
not  alter  the  system's  behavior. 

Ill)  The  Initial  probability  distribution  of  units  within  the 
compartments  Is  known  and  In  particular  the  cumulant 
generating  function  of  this  distribution  Is  denoted 

k(0j^,  02»  ••••  • 

1.5  Overview 

In  Cardenas  and  Matls  [1974],  the  stochastic  solution  to  the  n- 
compartment  catenary  and  mammillary  models  with  Irreversible,  time- 
dependent  transition  probabilities  Is  solved.  The  schematic  Illustrat- 
ing the  catenary  model  Is  shown  In  Figure  5. 


I 


I 

I 

I 


I 

1 

I 

I 

I 


Xoi(t) 


Figure  5 

n-Compartment , Irreversible  Catenary  System 


The  joint  cumulant  generating  function  of  the  stochastic  vector 
(X^(t),  X^Ct),  X^(t))  Is  found  In  terms  of  the  cumulant  generating 

function  for  the  Initial  distribution,  k(0j^,  9^)  • When  the 

Initial  number  of  units  Is  assumed  to  be  X^(0)  « X^,  1 •>  1,2,  ...,  n, 
the  form  of  this  generating  function  Is  shown  by  Hatls  [1970]  to 
correspond  to  a sum  of  multinomials  and  Poissons. 

In  a similar  manner,  the  Joint  cumulant  generating  function 
associated  with  the  mannlllary  model  Illustrated  In  Figure  6 Is  derived 
In  Cardenas  and  Hatls  [1974]. 


Figure  6 

n-Compartment,  Irreversible  Mammillary  System 


These  two  systems  are  combined  and  the  resulting  mixed  model  Is 
derived  In  Section  2 of  the  present  study.  Using  the  cumulant  generating 
function,  tlte  expectation  and  variance  for  each  conqpartment  as  well  as 
the  covr.rlance  between  any  two  compartments  are  found  as  a function  of 
time. 

This  model  is  generalized  further  In  Section  3 to  a coiq>letely 

general  Irreversible  model  In  which  a unit  In  compartment  n can  transfer 

St 

to  any  of  the  n-1  remaining  compartments  and  a unit  In  the  (n-1) 
compartment  can  transfer  to  any  of  the  n-2  other  compartments,  etc.  As 
In  Section  2,  the  cumulant  generating  function  and  the  first  and  second 


I 


I 

I 

I 


moments  of  the  stochastic  distribution  are  derived. 
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Section  4 contains  applications  to  reliability  analyses  using  a 
compartmental  approach.  Systems  of  components  In  various  redundant 
configurations  are  considered  and  standard  properties  such  as  the 
system  reliability,  the  mean  time  to  system  failure  and  the  variance  of 
the  llfelength  are  derived  as  functions  of  the  transition  rates. 

Other  areas  of  research  and  possible  applications  In  compartmental 
analysis  are  discussed  In  Section  5. 


I 


1 
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2.  THE  MIXED  MODEL 

In  this  section  a compartmental  model  is  considered  which  Is  a 
mixture  of  the  catenary  model  and  mammillary  model  discussed  In  Section 
1.5.  The  form  of  this  mixed  model  Is  Illustrated  In  Figure  7.  The 
stochastic  variable  Xj^(t)  Indicates  the  nisnber  of  units  in  the  1^^  com- 
partment at  time  t.  The  transition  rates  are  the  intensity  coef- 

ficients for  a unit  moving  from  compartment  j to  compartment  1 at 
time  t. 

The  exact  multivariate  distribution  for  the  stochastic  vector 
(Xj^(t),  X2(t) , ...,  X^(t))  will  be  determined.  The  method  employed  In 
this  derivation  was  first  used  for  time-dependent  compartmental  models 
by  Cardenas  and  Matls  [1974]  In  solving  the  n-compartment  catenary  and 
mammillary  models  previously  discussed. 

The  technique  consists  of  first  finding  a partial  differential 
equation  for  a generating  function  of  the  distribution.  This  partial 
differential  equation  Is  then  solved  and  the  multivariate  distribution 
Is  Identified  by  appealing  to  the  unique  form  of  the  generating 
function. 

One  method  for  finding  the  desired  partial  differential  equation  Is 
to  find  the  Kolmogorov  forward  equation  associated  with  this  model's 
single  event  probabilities.  From  the  forward  equation  the  partial  dif- 
ferential equation  follows  easily.  This  method  Is  Intuitively  appeal- 
ing. Essentially,  one  starts  %rlth  the  definition  of  a pattlal  deriva- 
tive and  proceeds  to  the  partial  differential  equation.  However,  this 
.ipproach  Is  tedious  and  a much  more  elegant  technique  was  devised  by 
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Bartlett  [1949],  which  allows  one  to  write  down  the  partial  differential 
equation  with  relative  ease.  A description  of  this  "random  variable" 
technique  can  also  be  found  In  Bailey  [1964]. 

2.1  Three-Compartment  Mixed  Model 

To  Illustrate  the  solution  to  the  mixed  model,  the  three-compart- 
ment model  depicted  In  Figure  8 will  first  be  solved  In  detail.  Then 
the  solution  for  the  n-compartment  mixed  model  will  be  sketched 
following  analogous  steps. 


^02^^^  ^10^^^  ^01^^^ 


Figure  8 

Three-Compartment  Mixed  Cstenery-Mammlllary  Model 
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2.1,1  The  Partial  Differential  Equation  for  the  Cuaulant  Generating 
Function 

All  of  the  possible  transition  probabilities  can  be  suanarlzed 
as  follows: 

Prob  {a  single  unit  enters  conpartnent  1 from 

outside  the  system  In  the  Interval  (t,  t+At)}- 
X^Q(t)At  + o(At)  for  1 - 1,2,3. 

Prob  {a  single  unit  moves  from  compartment  j 

to  compartment  1 in  the  Interval  (t,  t+At)> 
Xj(t)X^j(t)At  + o(At)  for  j - 1,2,3,  1 - 0,1, 2, 3. 

There  are  nine  events  with  probability  of  first  order  magnitude  of 
At.  All  other  events  can  be  Ignored  because  their  probability  Is 
o(At) . These  nine  events  Imply: 


Prob  {AX^(t)  • AX^Ct)  - k^,  AXjCt)  - kjl 

X^(t),  X2(t),  X3(t)}- 


X3(t)X^3(t)At  + o(At)  kj— 1,  k^M-1  1-1,2. 

X.(t)Xo .(t)At  + o(At)  k.--l,  J-1,2,3,  all  other  k*s 
^ ^ are  zero. 

X2(t)Xi2(t)At  + o(At)  k^"-!.  k3«+l  all  other  k's 

are  zero. 

X,Q(t)At  + o(At)  k,-+l,  j-1,2,3,  all  other  k's  are 
^ ^ ^ zero. 

where  AX^(t)  - X^(t+At)  - X^(t)  for  all  1. 
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Now  applying  Bartlett's  "randoa  variable"  technique  one  obtains  the 
partial  differential  aquation 


aMCe^.e^.e^.t)  2 -e,+e 

~±li~  ‘'"±3 


-e,+e  3M(e  0 0 t) 

’^>^13<‘> — hi  - 


3 -e  3M(0,,0,,0-,t) 

+ Z (e  ^-1)X  (t)  

J-1 


-0,+«,  3M(0,,0  ,0  t) 

+ (e  ^-l)X,,(t)  -i  ,1  3 


3 0 


+^Z^(e  ^-l)XjQ(t)M(0^,02,03,t) 


(1) 


for  the  Joint  moment  generating  function  of  Xj^(t),  X2(t)  and  X^(t) . 
After  dividing  (1)  by  M(0j^, 02*63,1)  and  rearranging  terms,  the  differ- 
ential equation  for  the  cuoulant  generating  function  K(0j^,02,03,t)  Is 
found  to  be 


3K(0^,02,03,t)  -0^ 


at 


- (e  -l)XQj(t) 


3K(0^.02*®3*t) 


-02-^1  1 
-1)Aj^2^‘>J 


3K(0^,02,03,t) 


30, 
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r -®3  2 -6+6  I 3K(6  e,,e,.t) 

-l)X„3(t)+  *-l)X3,(t>J L-f-^ 

2.1.2  The  Joint  Cu«ilant  Generating  Function 

Equation  (2)  la  a quasl-llnear  partial  differential  equation  which 
can  be  solved  using  characteristic  theory  (see  Ford  [1955]  or 
Garabedlan  [1964]).  The  characteristic  equations  associated  with  (2) 
are  determined  by  the  system  of  ordinary  differential  equations 


dt  °”1 ^^2 

1 " "®1  " ‘®2 

(1-e  ^)Xo.(t)  (1-e  ^)X_,(t)+(l-e  ^ bx,,(t) 


■®3  2 -®3'^l 

(1-e  ^)X,„(t)+  S (1-e  ^X,,(t) 


<K(e^,e2.e3,t) 

^ ®1 

Z (a  2-l)x,„(t) 


t 
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Letting  - (e  -1),  and  hence  - e de^^,  these  subsidiary 


equations  can  be  rewritten  as 


dVj  - V^XQj(t)dt, 

dV2  - [(V2-Vj^)^12^*'^  + V2Xo2(t)]dt, 

dV3  - ((V3-V2)X23(t)  + (V3-V^) X^3(t)  + V3XQ3(t) ]dt , 


t)  - V"]  «• 


This  system  of  differential  equations  can  be  solved  sequentially, 
obtaining. 


Xoi<z)dz 


V3  - C3  ezp< 


1 — JL 

T 

i 
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['*■  r *•  t 

/^Xi3(t2)exp|/^2Xoi(*)(ltf/^  (X23(a)+Xi3(*)+Xo3(*)  Id^dt^ 


J 


and 


K(03^,e2.03.t)  - C^-t«3/^^30(t3)exp|j^^[Xl3(z)+X23(*)+XQ3(z)]d*|, 


t.t. 


“ J / ^30^^3^^23^^2 


0 0 


)exp[/%Xu, 


(z)+XQ2(*)ld* 


f 3 ^ 


izjdtadt3 


■Vo^^20^V^12^V**p(/q\i(*>‘»* 

t t 

"Vo  J*** 


. ' • .'•.4  ■, 
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2“^3 


t»  t- 

+/  [Aj^2(’'^'^^02^*^  ^I^T,(z)+^7-»(z) 


'13'^' '"23' 


'j'itldtz’itj]. 


(5) 


Now  again  working  sequentially  with  (5)  and  substituting  back  the 
values  of  V^,  ^3  solve  for  the  arbitrary  constants 
^1*  ^2*  ^3  ^4*  Thi®  results  In 


Cl  - Uj^(0j,02*®3»^)  “ -Dp^i^t)  ’ 

®1  ®2 

^2  “ “2^®1’®2*®3’^^  “ -l)Pj^2(^^‘^(*  -l)P22(t)  > 

6 6 0 

C3  - 03(0^,02.03,0  - (e  ^-l)pj3(t)+(e  ^-l)p23(t)+(e  ^-l)p33(t)  » 


C4  ■ u^O^, 02,63,0  - K(0,,0,,0,,t)  - (e  -1)6,  (0 


3’ 


6 6 

- (e  ^-1)62(0  - (e  ^-1)63(0  » (6) 


where 


Pll(t)  - exp^-j^XQj^(s)dz^  , 


^/^Xo^(z)dzj 

P22<0  - exp|^/^[X32(*)+^02^*^^d*^  ^ 


»y 
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p^jCt)  = exp|^J^[Xj^^(2)+X22(z)+XQ2(z)]dzJ  , 

PiaCt)  » / XQj^(z)d^dtj^  T 

t t /'  t- 

Pl3(t)  » / / Xj^2(tj^)^23^*^2^®*PH 
0 t2  ® 

^2  *^1 

+ / X^-(t,)exp|-j’  ^[X  _(z)+X  (z)+X  (z)]dz 


'13 '"2 


13'  ' ''23'  ' 03' 


-/  X (z)d^dt2t 

t c t 

P23<t)  = / X23(t2>exp|-/  [X^^(z)+X23(z)+Xq^(z) ]dz 

h 

t r t > t t 

6j^(t)  ■ /^Xj^Q(t3)exp|-/  [XQ^(z)di^dt3+/^/^  ^20^^3^ ^12^^!^ 

expi-J  ^01^*^**^^1^^3 

^ ^3  H ^ 

t t ft 

+/  / X3Q(t3)Xj^3(t2>expM  ^[Xj^3<z)+X23(z)+Xq3(z)  ]dz 
0 t3  t3 

t -s  t t t 

to  ^ ^ 
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and 


r t t 

L 

-/  XQj^(z)d^dtj^dt2dt^  ? 

62(0  =•  J X2Q(t2)exp^/  [Xj^2(’')+^02^^^^'*^‘^*^3 

t t /-  t„ 

+/  J X^(t^)X23(t2)exp|-J  [Xj^^(z)+X22(z)+Xq2(z)  ]dz 

0 t3  K. 

-j  [Xj^2^*)'*'^02^’'^ 

•^2 

«3(t)  - / X^Q(t2)exp|^/  [Xj^^(z)+X22(z)+XQj(z)]d^dt2. 


These  four  arbitrary  constants  (6)  are  the  Independent  Integrals 
which,  according  to  characteristic  theory,  can  be  related  by  the 
functional  relationship 


u^(0j^, 02* ®3» ’ **2^®1’ ^2* ^3* ’ 


Uj(0j^, 62 , Oj, t) ) ' 


(7) 


The  exact  functional  fora  of  ^ can  be  Identified  If  the  Initial 
conditions  of  the  aodel  are  known.  For  this  purpose,  assuote  X^(0) , 
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X2(0)  and  X^CO)  have  the  joint  cumulant  generating  function  k(0j^,02.02) f 
l.e.  K(0^, 02,63,0)  - k(0j^, 02. 63). 

Also  when  t * 0 the  system  of  equations  (6)  becomes 


®2 

“ e -If 

63 

'^3^®1’®2’®3’°^  = e -1  ) 


u^(6i, 62,03.0)  = k(0^,02,03). 


(8) 


So  when  t * 0,  eqxxation  (7)  reduces  to 


h ®2  ®3 

k(0i,02,63)  = i(<(e  -1,  e -1,  e M) 


(9) 


Now 


letting  y^  ■ e -1  for  i ■ 1,2,3,  (9)  can  be  rewritten 


k(in(yj^+l),^n(y2+l),<.n(y3+l))  - <Kyj^,y2.y3) 


(10) 


Equation  (10)  gives  the  functional  form  of  v In  terms  of  the  Initial 
Joint  cumulant  generating  function.  Using  this  form  for  ij>  , equation  (7) 


becomes 
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«,n(u2(6j^,02,e3,t)+l) , 

An(u2(ej^,02.03,t)+1))  . (11) 

The  solution  to  (2)  is  now  obtained  by  substituting  (6)  Into  (11) 
and  solving  for  K(0j^,02,02»t)  • This  results  In 

r 20, 

K(0^, 02,03, t)  - kjln(l+(e  -1  )pj^j^(t)),  ln(l+  E (e  -l)p^2(‘'>^« 

3 0 -1 

ln(l+  E (e  ^-l)p  (t)) 

i-'l  J 


3 0 

E (e  ^-l)6,(t) 

i-1  ^ (12) 


where  the  and  6^(t)  are  defined  in  (6). 


2.1.3  The  Moments  of  the  Mixed  Three-Compartment  Distribution 

The  moments  of  the  distribution  can  now  be  found  by  differentiating 
the  cumulant  generating  function  (12)  with  respect  to  the  appropriate 
parameters  and  evaluating  the  partial  derivatives  at  0^^  ■ 0 for  all  i. 
DefindLng  and  to  be  the  mean  and  variance  of  Xj^(O)  for  1-1,2 
and  3 and  to  be  the  covariance  of  Xj^(O)  and  Xj(0),  1,J  - 1,2,3, 
i J,  then  the  moments  are  given  by 
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E(X^(t)) 

V(X^(t)) 

cov[X^(t),  X^(t)] 


3 

£ 

k.-l 


3 3 

k=«l  4-1  ^ 


3 3 

k*i  Jl-j  ^ 


^ W,P4,(t)P,«(t) . 


X,-max[i,  j ] 


(13) 


These  results  are  a generalization  of  the  2- compartment  catenary 
model  and  the  3-compartment  mammillary  model  of  Cardenas  and  Matls 
[1974].  See  Appendix  B for  the  details  of  this  derivation. 

2.2  The  n-Compartment  Mixed  Catenary-Mammillary  Model 

2.2.1  The  Partial  Differential  Equation  for  the  Cumulant  Generating 
Function 


The  3-compartment  mixed  model  will  now  be  extended  to  n compart- 
ments. The  derivations  of  the  joint  cumulant  generating  function  and 
the  moments  of  the  distribution  are  completely  analogous  to  those  of 
the  3-compartment  case.  Figure  7,  page  17,  Illustrates  the  n-compart- 
ment  mixed  model  under  consideration. 
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The  transition  probabilities  are: 

Prob  {a  single  unit  enters  compartment  1 from 
outside  the  system  in  the  Interval 
(t,t+At)}  - X^^(t)At  + o(At)  for  1-1,2, ... ,n, 

Prob  {a  single  unit  moves  from  compartment  j to 
compartment  1 In  the  Interval  (t,t+At)}  - 
Xj(t)Xij(t)At  + o(At)  for  j=l,2,...,n 
and  1-0, 1,2 n. 

There  are  4n-3  events  with  probability  of  first  order  magnitude 
of  At  and  the  probability  of  these  events  will  be 

Prob  (AXj^(t)-kj^,  AX^CO-k^ AX^(t)-k^|Xj^(t) 

Xn(t)X^n^t)At-K)(At)  k^-fl,  k^— 1 1-1,2, ...  ,n-l 

all  other  k's  are  zero. 

Xj^(t)Xoi(t)At+o(At)  k^— 1 1-1, 2,..., n 

all  other  k's  are  zero. 

Xi(t)Xi_i  ^(t)At+o(At)  k^— 1 l-2,3,...,n-l 

all  other  k's  are  zero. 

Xio(t)At+o(At)  k^-+l  1-1,2, ... ,n 

all  other  k's  are  zero. 

The  differential  equation  for  the  Joint  cumulant  generating  function 
can  now  be  written 
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3K(i.t) 

at 


- (e 

Kin  (t)i^^^^ 

n-1 
+ Z 
1-2 

r •®i'*^i-i 

l(. 

r - 

•6  n-1  6 

+ (e 

”-l)X  (t)+  Z (e 

0“  1-1 

J 


n e 

+ Z (e  M)A.-(t), 
1-1 


(14) 


where  K(6,t)  denotes  KCOj^.e^, . . . ,0^,t)  . 

2.2.2  The  Joint  Cuaulant  Generating  Function 

The  subsidiary  equations  associated  with  equation  (14)  are 


1 


d9. 


d0. 


-e. 


-e. 


/ 6 ••0  — 

(l-e  (1-e  Sx^_j^^^(t)+(l-e  "^)X(,^(t) 


do 


■ n 1-2,3. . . . ,n-l. 

I (e  ^-l)X,„(t) 

1-1 


(15) 


Letting  - c -1  for  1 - l,2,...,n  the  subsidiary  equations  can  be 


rewritten 
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■IVl  - VjXjj(t)dt  , 


'■''l  ■ ‘-2.3 n-1, 


[»nN)a<«+"yVW<')]' 


dK(0,t)  - Z V.X..(t)dt. 
1-1  ^ 


This  system  can  now  be  solved  sequentially  for  V , 1-1, 2,..., n 
and  K(e,t)  In  terms  of  the  transition  rates  and  irt-1  constants  of 

Integration.  Solving  this  system  for  the  constants  of  Integration  and 


replacing  by  (e  -1)  one  obtains 


Cl  - Uj^^.t)  - (e  -l)e3qp{-gj^(t,0)}  , 


9,  t 


1 

+ (e  -1)/ 


04 

Cl  - “id.t)  - (e  M)exp{-gj^(t.O)} 


04_i  t 

. mv  / 


1-2  e,  t t t 


^1-1^*  "^^Vt  *’*^t  

° *=1-1  ^1+1 


33 


4 


1-1  1-1 

exp(-g,(t,t,)-  ^ «k(‘Tc-rV-8i<Vr°>)  "/‘tc 
K"*+l  k*l 

for  1 - 3,4,...,n-l, 


-l)exp{-g^(t,0)} 


n-2  6,  r.t 


^ m=l  Vt”*  ^t  ^^■Ha.n^Vl^^t.l-hn^*'l ‘'l-Hn-1^ 


n-1  "1+1 

1+m 


*Tm  £+01-1  -1 

«p(-8j(t.t,)-J^^g^(t^.j.t^)-8„(t„.j.0»  -wJ  , 

t t 

n 

t 1 t 

£-2  0 tj^j^^J+£,0^*^n^*^J,J+£^^J’*“^J+£-l^**P^"*j^^»^j^ 

j+£-l  j+l-1 

"k-j+l*k^*^‘‘-l’'^^"*J+*^’^J+Jl-l’*^n^^  klj 

■"V,  *nO<',>S"‘'»-l’"''‘'*j‘"’Vl>-*.<'i.-l'V><‘Vl'“n 


m 
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n j+1 

J+A  j+A-1 

«I.(-8j(t.tj)-  Z ^(t^.i,t^)-^(t^_l,tn)l  n dVt„.jdt„, 


k-J+1 


k-3 


(17) 


A 

where  gj^(x.y)  - / Xq^(z)<1z  , 


gj^(x.y)  - / t^oi^^^'‘'^l-l  ’ t-2,3,...,n-l  , 

y 


X n-1 

g^(x,y)  - / S X.  (z)dz 


n 


y J-o 


Jn' 


and 


These  n+1  constants  can  now  be  related  by  some  function  ^ such  that 


- ♦(Uj^(0,t) u^Ci.t)) 


(18) 


where  ♦ is  determined  by  the  Initial  conditions.  Assume  the  Initial 
distribution  Xj^(O),  X2(0)  , . . .X^(O)  has  the  joint  cumulant  generating 


function 


I 

I 

I 

I 


K(02’®2*  " ’ * " k(0j^, 02 » . « . » 0^) 


(19) 
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Now  one  can  verify  that  the  joint  cvimulant  generating  function  of  the 
stochastic  vector  Xj^(t)  , . . . ,X^(t)  at  any  time  t Is  given  by 


K(01.02 ®n**=^  “ k^An[l+(e  M)p^^(t)l,  to[l-f^Z^(e  J • 


J ®1 

. to[l+  E (e  -l)p,,(t)l 


1-1 
n 0 


Ij 


no  > n 0 

ln[l+  (eM)p  (t)]^+  Z (e  J-l)6,(t), 

1-1  ^ J J-1  J 


(20) 


where 


Pil(t)  =■  exp{-g^(t,0)}  I 


t t 


1 = 1, . . . ,n  > 


L/  "•/  rj^^(t^,...tj^_j^)exp{-gj^(t,tj^) 


■1-1  ‘m 


1-1 


1-1 


~ ^ 1’ ^ * 
k-m  *^  *^  ^ *^  1-  k-i 


for  1 - 2,3,...,n-l;  I - 1,2,. ..,1-1. 
t 

r > 

0 

n-1-1  t t t 


Ptn<">  “ /o\n<Vl>®*P^-«l<‘=*Vl^-«n<Vl*0^^‘*>^n-l 

^ m-1  Vt  ^t-Hn,n^Vl^^l,l-Hn^*^l’---*^l+iB-l) 

n-1  't+1 

t+B  l+n— 1 


for  I - 1,2, . . .n-1. 


■irr-«r- 


fcftln  k ■ 
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and 


t t 


expl-gj(t.tj)-gj^j(tj,t„) )dtjdt^ 

U-J-1  t t t 

^ a-2  Vt  ”*^t  Va-1) 

n J+1 

j+H-1  j+l-1 

'no ‘V  ' Jn' 'n-l> < ' ■ Vl> -*n< Vl • 'n>  >'*'n-l'‘'n 


* g-1  Vt  "'^t  V«-l*  Vx.n^Vl' 

n J+1 

j+H  J+A-1 

d*p{-gj(t,tj>-  £ ek(Vi.V-8„(t„.i.t„)l  It  ■‘Widt,^ 


k^J+1 


k-j 


where  gj^(x.y)  and  (z^^, . . . ,Zj_  j,)  are  defined  in  (17). 

2.2.3  The  Moments  of  the  Mixed  n-Compartment  Distribution 

The  moments  derived  from  this  cumulant  generating  function  (20)  are 
the  logical  extension  of  those  derived  for  the  3-compartment  model.  As 
before,  define  and  to  be  the  mean  and  variance  of  X^(0)  for  all 
1,  and  0^^  to  be  the  covariance  of  X^(0)  and  Xj(0)  for  all  1 and  J, 

1 J.  (See  Appendix  B for  Che  details.)  One  obtains 


I 
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E(X^(t)) 

V(X^(t)) 


^ ^k^ik 
k-i 


n ^ n 
E 7 E a 
k-lLi-l 


(t)-W^(t),  i-1,2 


«^(t), 


i-1,2 


f • • • • 


n 


n n n 

cov(X  (t),X  (t))  -EE  ^ * 

^ J k-1  i-j  i-max[l,j]  ^ 

i,j-l,2,...,n  i^j,  (21) 

where  P^j(t)  and  6^(t)  are  defined  in  (20)  for  all  1 and  j. 

Higher  order  moments  may  be  obtained  by  a straightforward 
extension  of  this  procedure  to  higher  derivatives. 


2.3  A Special  Case 


It  is  frequently  true  that  the  initial  number  of  units  in  each 

compartment  is  known.  For  exeimple,  in  the  gastrointestinal  model  cited 

in  Section  1.2,  it  was  known  that  the  rumen  contained  4,000  beads  at 

time  zero  while  the  two  other  compartments  contained  no  beads. 

If  the  initial  distribution  is  known  to  be  X^^(O)  - for 

1-1, 2,..., n then  K(0,O)  - k(0)  - X, 0,+X_0,+. . .+X  0_.  Now  the  cumulant 
— — xizz  nn 

generating  function  can  be  written  from  equation  (20) . This  becomes 


n J 

K(0,t)  - E X.inll+  E (e  -l)p.,(t)  1 
J-1  J i-1 


+ E (e  ^-1)0, (t) 

J-1  ^ 
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where  the  P^j(t)  and  6j(t)  are  defined  In  (20). 

Metis  [1970]  discusses  a cuaulant  generating  function  of  this  form 

extensively  and  shows  It  Is  associated  with  a sum  of  n Independent  multi- 
nomial distributions  and  a sum  of  n Independent  Poisson  distributions. 

Establishing  this,  one  can  then  appeal  to  the  unique  form  of  the  gener- 
ating function  and  realize  a physical  Interpretation  of  the  and 

6j(t)  parameters.  This  Is 

p^j(t)  “ Prob  {a  unit  which  was  In  compartment  j at  time 
t"0.  Is  In  compartment  1 at  time  t}, 

6j(t)  =«  The  expected  number  of  units  In  compartment  J at  time  t 

which  were  not  In  the  system  at  time  t * 0. 

This  result  Is  very  important  In  the  later  development  of 
applications  using  compartmental  analysis. 
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I 
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I 

I 

I 

I 
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3.  THE  GENERAL  IRREVERSIBLE  TIME-DEPENDENT  MOm 

The  n-compartment  mixed  catenary-mamillary  model  derived  In  the 
previous  section  Is  further  generalized  In  this  section.  The  assumption 
Is  made  that  a unit  In  the  first  compartment  can  transfer  to  any  one  of 
the  remaining  n-1  compartments  and  a unit  In  the  second  compartment  can 
transfer  to  any  one  of  the  remaining  n-2  compartments  and  so  on.  In 
addition,  a unit  can  enter  or  leave  the  system  through  any  compartment. 
This  arrangement  Is  Illustrated  In  Figure  9.  Under  the  assumptions  In 
Section  l.A,  such  a structure  Is  the  most  general  possible  Irreversible 
time-dependent  model. 


3.1  The  n-Compartment  General  Model 

3.1.1  The  Partial  Differential  Equation  for  the  Cumulant  Generating 
Function 

The  possible  transition  probabilities  can  be  summarized: 

Prob  {single  unit  moves  from  compartment  j to 
compartment  1 In  the  Interval 
(t,t+At)}  - Xj(t)X^j(t)At+o(At) 

for  J«2,3,...,n  and  1-1,2 j-1. 


41 


J 


t 


Prob  {single  unit  In  compartment  J leaves  the  system 
In  the  Interval  (t,t+At)}  ■ 

Xj(t) AQj(t)At+o(At)  for  all  J. 

Prob  {single  unit  enters  compartment  j from  outside 
the  system  In  the  interval  (t,t+At)}  - 

XjQ(t) At+o(At)  for  all  j. 

2 

In  this  model  there  are  (n  +3n)/2  events  with  probability  of  first 
order  magnitude.  These  events  lead  to  the  following  probability 
statement. 


Prob  {AX^(t)*k^,AX2(t)=k2, 
Xj(t)Xoj(t)At+o(At) 

« Xj(t)X^j(t)At+o(At) 

XjQ(t)At+o(At) 


. . .AX^(t)-k^|x. (t) ,X,(t) X^(t) } « 

n n 1 / n 

kj«-l  for  j«l,2,...,n 

and  all  other  k's  are  zero. 

kj«-l,  k^»+*l  for  j«2,3,...,n  and 

1-1,2,. ...j-l. 

kj»H-l  for  j*l,2,...,n 

and  all  other  k's  are  zero. 


Now  proceeding  as  in  the  mixed  model  the  multivariate  distribution 

of  the  n-vector  (X, (t),  X.(t),...,X  (t))  will  be  determined  by  finding 
J.  z n 

its  cumulant  generating  function  K(0j,02, . . . ,9^,t) . Denote  this  by 
K(^, t).  The  partial  differential  equation  for  K(^,t)  can  be  shown 
to  be 
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1 J-2 


J-1  -e  +e 

E (e  J M)X..(t) 
1-1 


(22) 


3.1.2  The  Joint  Cumulant  Generating  Function 

The  details  of  the  solution  to  (22)  are  found  In  Appendxx  A.  The 
method  is  the  same  as  that  used  In  solving  (2) . The  solution  is 


K(0,t) 


r ®1  ^ \ 

= k|ln[l+(e  M)p^^(t)],  ln[l+  E (e  1 » ' • • « 

J 01  n 0 - 

ln[l+  E (e  -l)p.,(t)] ln[l+  E (e  M)p,  (t)  ] 

1-1  ^ 1-1  ^ 


n 0 


+ E (e  ^-1)6  (t) 
.1=1  ^ 


(23) 


where  K(6,  ,e«, . . .,0  ,0)  - k.(6.,  ,9_, . . . ,0  ), 
1 4 n i.  4 n 


Pll(t)  - exp{-h^(t,0)},  1-1,2,. ..,n. 


?1  i+i(t)  “ / ^^j^(t^)exp{-h^(t,t^)-h^^j^(t^,0)}dt^,  1-1, 2, . . . ,n-l. 


Pij(t)  = J Iij(tj_j^)exp{-h^(t,tj_j^)-hj(tj_j^,0))dtj_j^ 


j-1  t t 

1-1+1  0 tj_j^  1*.  J J.  1 «.  i 
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^ t t t 

^ ^ ^ ^ ’"^  -1^ 
l<tj«2<--<WJ  » Vl  '*2-1  ^ ' 


m-2 

n (X  (t  ))x 

i-l  1’  1+1  1+1  m-l’ 


j(tj_^)exp{-h^(t,t^^_^) 


m-2  m-1 

— ^ h.  (tj  i»^j  -1  »o) } n 

1-1  *■!  i"^  i+r^  * 1-1  *'l"^ 

t t t 

"'K  ^ij^‘^l’*^l+l’"'’  tj_j^)exp{-h^(t,t^) 

° ‘^J-1  *^l+.l 

j-1  J-1 

" ^ h,(t  .,t  )-h.(t  0)}  n dt,p 

1-1+1  *•  *•  1 J-X  *■ 

j-3,4,...,n  and  1-1,..., j-2. 


t t t 

6 ,(t)  » / X , ^(t  )exp{-h  ,(t,t  ) }dt  +/  / X _(t  )X  , (t  .) 

n-1  n-1,0  n n-1  * n n ^ nO  n n-l,n'  n-1 

n 

exp{-h  ,(t,t  J-h  (t  ,,t  )}dt  ,dt  . 
n-x  n— 1 n n-l  n n-1  n 


V2<'>  • ^gV2,0<V“>'<-''«-2"-'n>>‘''n 

t t 

■*■/  / 1 9 « 1^^..  -)exp{-h  ,(t,t  J 

^ n-l,u  n n-2,n-l  n-2  n-2  n-2 

n 

t t 

-h  ,(t  -,t  )}dt  ,dt  +f  f X -(t  )X  , ^(t  ,) 

n-1  n-2  n n-2  n nO  n n-2,n  n-1 


exp{-h  ,(t,t  )-h  (t  ,,t  )}dt  ,dt 

n-2  n-1  n n— i n n-1  n 


t t t 


+/  / / ^ n(t  ) A , , (t  .)  X , (t  J 

' nO  n n-2,n-l  n-2  n-l,n  n-1 


0 ^ ^ y »•  A ft*  ftft  y U ftft 

n n-1 

exp{-h  ~(t,t  ,)-h  ,(t  ,t  ,)-h  (t  ,t  )}dt  _dt  ,dt  . 

n-z  n-Z  n-i  n-Z  n-i  n n-l  n n-z  n-i  n 


n t t 


U l*j‘rX  U t 


exp{-hj(t,t^_^)-h^(t^_^,t^)}dt^_^dt^ 
n-j  j+A-1  t t t 

_ t t t t 

■*■  / / / / **"^  I -1^ 

J<iT<A-<.V.<*,  ,<n  0 t t.  , t.  ^*1 

12  m-1  n l^-l 


n A«  « (tj  )exp{-h,  (t,t  £ h (t  _i»^«  _i^ 

i-1  ^l’*'l+l  ^1+1  ^ J *•!_  J-  \ \ ^ i 


m-1  m-1  1"1  1 

"'K  L ^nO^V^jn^^J’*^J+l Vl^ 

” n J+2  J+1 

exp{-hj(t,tj)-  £ n dtj^. 
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where 


X i-1 

h (x.y)  - / Z X..(z)dz 
y J-0 


and 


This  result  is  a generalization  of  the  n-compartment  mixed  model 
derived  in  Section  2.2. 


3.1.3  The  Structure  of  p^j(t)  and  6j(t) 


It  is  evident  p^j(t)  and  6^(t)  are  lengthy,  complex  expressions 
and  this  complexity’  Increases  rapidly  with  an  Increase  in  the  number  of 
compartments. 

But  the  structure  of  these  terms  can  readily  be  seen  by  considering 
their  relationship  to  the  binomial  coefficients. 

As  pointed  out  in  Section  2.3,  p^j(t)  is  Che  probability  a unit  in 
compartment  J at  time  zero  will  be  in  compartment  i at  time  t.  Since 
all  Che  possible  paths  from  compartment  j to  compartment  1 are  mutually 
exclusive  events,  p^j(t)  must  be  the  sum  of  the  probability 
contributions  from  each  of  these.  The  probability  of  each  path  is 
expressed  as  a multiple  integral  where  the  multiplicity  is  equal  to  Che 
number  of  transfers  along  that  path.  A given  transfer  is  either  a part 
of  a path  from  compartment  J to  compartment  i or  it  is  not,  and  it  is 
this  binomial  aspect  which  gives  rise  to  the  binomial  coefficients  of 
Figure  10.  This  figure  is  Che  well  known  arrangement  of  the  binomial 
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coefficients  known  as  Pascal's  triangle.  The  Integral  structure  of 
p^j(t)  is  given  by  the  rows  of  the  triangle.  The  diagonals  counting 
from  the  left  correspond  to  the  multiplicity  of  the  integral.  For 
example,  if  J-i-k,  then  the  k*^**  row  indicates  will  have  the  form; 

^one  single  integral^  + ^ (k-1)  double  integrals^  + . . . 

^{(^0  integrals^  + ...  + ^one  k-fold  integral^. 

k-1 

This  is  a total  of  2 Integrals  to  evaluate. 


12  1 
13  3 1 


single  Integrals 
— double  integrals 
triple  Integrals 
four-fold  integrals 


Total  no. 
of  terms 


1 4 6 4 1 


m-fold  Integrals 
k-fold  Integrals 


1 k-1 


a- 


(n-l)-fold  Integral 


n-1  1 n-2 


Figure  10 

Structure  of  uaing  Pascal's  Triangle 
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In  an  application  of  this  model  to  reliability  in  Section  4 the 
calculation  of 


E P.„(t) 

i-1 


is  required.  This  could  contain  as  many  as  2°  ^-1  Integrals  to 
evaluate  I 


-f-". 


The  structure  of  determined  in  a similar  way  using 

Pascal's  triangle.  For  a model  containing  n compartments  the  expression 
for  contains  2°  ^ terms.  The  number  of  integrals  of  a given  mul- 

tiplicity can  be  read  from  Figure  11.  For  ^j(t)  the  (n-J+1)^^  row  of 


Row 

1 

2 

3 

4 

5 


^ — single  Integrals 
1 double  integrals 

1 1 ^ triple  Integrals 

12  1 four-fold  Integrals 

13  3 1 

14641  m-fold  integrals 


n-J+l  1 n-J  •••  •••  1 


n 1 


n-1' 
m— 1| 


(n-J+l)-fold  integrals 


n-fold  integral 


6^(t) 

'n-2<^> 

«n-3<^> 

V4^^> 


«j(t) 


6^(t) 


Figure  11 

Structure  of  6^(t)  using  Pascal's  Triangle 


Pascal's  triangle  indicates  the  expression  will  have  the  form; 


1 
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^one  single  integral^  + ^(n-J)  double  Integrals^ -t-  ... 

+ m-fold  Integrals^  + ...  + ^one  (n-j+l)-fold  integral^ 


3.1.4  The  Moments  of  the  General  n-Coapartaent  Distribution 


Define  and  to  be  the  mean  and  variance  of  X^(0)  and  to 
be  the  covariance  of  X^(0)  and  Xj(0).  Then  the  moments  of  the  general 
n-compartment  distribution  can  be  written  directly  from  (21) . Observing 
that  the  cumulant  generating  functions  (20)  and  (23)  are  Identical  In 
form,  the  moments  will  also  be  Identical  In  form.  Hence 


E(X^(t)) 

V(X^(t)) 


k»i 


n r n 
I 0 
k-iVl-i 


(t)  + 6j^(t)  1-1,2 n 


1-1, 2,..., n 

cov(X  (t),X  (t))  - E / Wj)l(t)p  (t)-  E 

k-1  1-j  ■'  l-max[l,j]  ■' 

l,j-l,...n  ii*j  (24) 


where  Pj^j(t)  and  d^(t)  are  defined  in  (23)  for  all  i and  J. 

Higher  momenta  may  be  obtained  by  a simple  extension  of  this 
procedure  to  higher  derivatives.  See  Appendix  B for  the  details  of 


this  derivation. 
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3.2  The  Four-CompartoenC  General  Irreversible  Model 


In  this  section,  the  formulae  for  the  general  n-compartment  model 
are  shown  for  the  special  case  when  n>4.  The  cumulant  generating 
function  becomes 


“0  2 0 

K(0j^,92,93,e^,t)  - k to[l+(e  ^-l)pj^^(t)  ] , ta[l+  E (e  ^ ♦ 

3 9 4 9 1 

toIl+^E^(e  -l)p^3(t)],  ln[l+^Z^(e  -l)p^^(t)]J 

4 9 

+ Z (e  ^-1)6  (t)  (25: 

j-1  J 


where 


Pll(t)  - exp{-h^(t,0) } 1 « 1,2, 3,4, 


Pl2<t)  - / 

t t t 

Pl3(t)  -/  Xj^3(t2)exp{-hj^(t,t2)-h3(t2,0)  }dt2+J  / ^i3(tj^.t2) 

exp{-hj^(t,tj^)-h2(tj^,t2)-h3(t2,0)  }dtj^dt2, 

t t t 

Pl4(t)  “ /^^i4^t3)exp{-hj^(t,t3)-h^(t3,0)}dt3+J  / ^24^*^3^ 

exp{-hj^(t,tj^)-h2(tj^,t3)-h^(t3,0)  }dtj^dt3 


^I3^t2)^34(t3)exp{-hj^(t,t2)-h3(t2,t3)-h^(t3,0)} 

dt2dt3 
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t t t 


-h^(t3,0)  Jdtj^dt^dt^, 


P23(t) 


2^  ^ 2'  3^  2*  ' **  2* 


P24(t) 


t t t 

^Q^24^‘^3^®*P^"*^2^‘^**^3^"N^^3*°^^‘*‘^3  ^24^‘^2’*^3^ 

exp{-h2(t,t2)-h3(c2,t3)-h^(t3,0) }dt2dt3. 


P34<t) 


/ ^34(t3)exp{-h,(t,t,)-h^(t,,0)}dt,. 


The  6^(t)  are  given  by 


«l(t) 


3 t t 


exp{-hj^(t,tj^)-h^_^j^(tj^,t^)  }dtj,dt^ 


t t t 


"^/  / / ^3q(^4^  Tj^3( t j^, t2) cxp{— hj^( t , tj^)— h2(t 12) 


0 t2 


t t t 


-h3(t2,t^)}dt^dt2<it^ 


+//  / ^4Q(b^)  ^24  ^^3^  ®*p{“bj^(  t , t j^)  ~^2^^1*  ^3^ 

0 t / t ^ 


t t t 


t \o'VH3<'2>*34<'3> 


0 t,  tj 


exp{-hj^(t,t2)-h3(t2,t3)-h^(t3,t^))dt2dt3dt^ 
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' 4.  SOME  APPLICATIONS  OF  COMPARTMENTAL  MODELS  TO  RELIABILITY 

I 

4.1  Introduction 

The  examples  in  Section  1.2  Indicate  the  units  within  the  compart- 
ments can  be  Interpreted  In  many  ways.  Recall  these  were  Inulln  within 
a rabbit,  people  within  an  organization,  beads  within  a gastrointestinal 
tract  and  carcinogenic  particles  In  a manmlllary  system.  In  reliability 
theory,  a "unit"  may  be  any  part  of  a system.  It  need  not  be  basic  or 
undecomposable . It  refers  only  to  a component,  however  complex,  whose 
reliability  is  of  Interest  to  the  overall  reliability  of  the  system. 

If  a unit  begins  operation  at  t ~ 0 and  failure  occurs  at  t ■>  t^,  then 
t^  Is  said  to  be  the  lifetime  of  that  unit.  Suppose  t^  Is  a random 
variable  with  distribution  function 

F(t)  •>  Prob  {tjj<t}. 

The  reliability  of  the  unit  at  time  t Is  the  probability  of  failure  free 
operation  In  the  Interval  (0,t)  and  the  reliability  function  Is  defined 
to  be  R(t)  - l-F(t) . 

The  distribution  function  Is  usually  assumed  to  be  an  intrinsic 
property  of  a unit.  And  all  such  units  are  assumed  to  be  homogeneous 
In  the  sense  that  their  life  lengths  are  distributed  Identically.  For 
this  to  occur  the  conditions  under  which  "Identical"  units  are  used 
must  be  homogeneous  otherwise  their  reliabilities  would  differ  due  to 
extrinsic  variables. 


I 


53 


In  this  section,  the  environments  in  which  a set  of  units  function 
are  assumed  to  be  (possibly)  heterogeneous.  The  function  X(t)  is 


assumed  to  reflect  the  operating  environment  of  a unit  at  time  t,  and 
X(t)At  can  be  Interpreted  for  small  At  to  be  the  probability  that  a unit 
which  has  functioned  without  failure  up  to  the  Instant  t will  fail  in 
the  interval  (t,t+At).  More  precisely,  X(t)At+o(At)  Is  the  conditional 
probability  of  failure  at  time  t given  the  unit  has  operated  success- 
fully up  to  that  instant. 

Letting  R(t,t+At)  denote  the  probability  that  a unit  has  functioned 
without  fail  during  the  interval  (t,t+At),  it  follows  that 


R(t,t+At) 


R(t)  • 


Hence  the  probability  of  failure  during  (t,t+At)  is 


A(t)At+o(At)  * 1-R(t,t+At)  * 


Therefore, 


hence 


^ ■ 


R(t)  l-F(t) 


provided  F'(t)  exists  and  F(t)<l. 
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Solving  (26)  for  R(t)  one  obtains  the  reliability  fimctlon  In 
terms  of  X(t) ; 

c 

R(t,  = exp{--j  X(z)dz}.  (27) 

0 

The  function  X(t)  Is  of  fundamental  Importance  In  many  fields  and 
Is  known  by  a variety  of  names.  In  actuarial  work  it  Is  the  "force  of 
decrement"  or  "force  of  mortality"  and  the  demographer  refers  to  it  as 
the  "age-speclf Ic  death  rate" . In  vital  statistics  and  extreme  value 
theory  It  is  the  "intensity  function".  In  reliability  theory  X(t)  is 
the  hazard  or  failure  rate.  Barlow  and  Prochan  [1965]  discuss  the 
role  of  the  hazard  function  In  reliability  analyses.  Barlow  et  al. 
[1963]  derive  many  useful  properties  of  R(t)  assuming  only  a ironotone 
hazard  function.  A listing  of  system  properties  using  the  monotone 
assumption  Is  found  In  Prochan  [1966]. 

The  Interpretation  of  hazard  rate  as  a conditional  probability 
corresponds  exactly  to  that  of  transition  rate  discussed  In  Section  1, 
provided  "leaving  a compartment"  Is  Interpreted  as  failure.  Hence,  the 
transition  rates  regarded  as  hazard  rates  In  a reliability 

setting  and  the  failure  distribution  can  be  found  In  terms  of  the  hazard 
rates  using  equations  (26)  and  (27) . 

The  hazard  rate  can  sometimes  be  decomposed  Into  constituent  parts 
having  particular  physical  Interpretations.  An  example  of  this  Is  seen 
In  a competlng-rlsk  model  where  a unit  can  have  several  Independent 
modes  of  failure.  The  hazard  rate  for  the  unit  Is  the  sum  of  the 
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individual  hazards  associated  with  each  risk.  Such  a system  is 
considered  in  Section  4.3.  There  are  also  situations  in  which  the 
hazard  rate  cannot  be  conveniently  decomposed  but  instead  the  failure 
distribution  is  written  as  a composition  of  distribution  functions. 
This  is  called  the  mixed  distribution  model  and  the  overall  failure  is 
written  as  a convex  sum  of  failure  distributions.  Kao  [1959]  used 
this  device  to  model  the  combined  failure  due  to  catastrophic  failure 
and  wearout  failure  as  a convex  sum  of  two  Weibull  distributions. 

The  mixed  distribution  model  allows  one  to  have  distinctly 
different  hazard  rate  functions  over  different  time  periods  in  the 
life  of  a device.  However,  in  this  model  the  time  domain  is  still 
referenced  to  the  moment  the  unit  is  turned  on. 

Mann  et  al.  [1974]  discuss  the  general  distribution  of  time  to 
failure,  wherein  a unit  can  fall  due  to  random  failure  (time- indepen- 
dent) or  due  to  wearout  or  Intrinsic  failure  (time-dependent) . The 
hazard  rate  then  has  the  form 


X(t)  • X , + X,  ^ ^ , (t). 

random  intrinsic 


(28) 


In  Section  4.2,  the  hazard  functions  will  be  of  the  form 


X(t)*X  +X  , .(t). 

random  extrinsic 


Unlike  the  wearout  hazard  function,  the  extrinsic  or  environmental 
hazard  function  may  be  referenced  to  a time  other  than  the  turnon  time 
of  Che  unit.  The  random  hazard  may  be  due  to  intrinsic  or  extrinsic 
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forces  and  of  course  has  an  exponential  failure  associated  with  It. 

Davis  [1952]  Is  a frequently  used  reference  for  Justifying  an  assumption 
of  exponential  failure  for  a wide  variety  of  systems. 

Some  results  of  Sections  2 and  3 will  now  be  considered  in  a 
reliability  context.  To  place  these  results  In  p>ersi>ectlve,  It  is 
necessary  briefly  to  discuss  some  fundamental  svstem  structures  and 
concepts  encountered  in  system  reliability  analyses.  The  reliability, 
R(t) , of  a system  can  be  Increased  in  many  ways  but  these  can  be 
generally  classified  as  either  redundancy  or  repair.  The  Irreversible 
compartmental  model  can  be  thought  to  represent  a redundant  system  with 
backup  components  in  a cold  standby  mode.  Hot,  warm  and  cold  standbys 
are  differentiated  according  to  the  manner  in  vdilch  the  standby  is 
loaded.  Hot  standbys  are  loaded  precisely  the  same  as  the  operating 
unit,  hence  are  assumed  to  fail  by  the  same  failure  law  as  the  opera- 
tional unit  even  though  not  in  use.  A warm  standby  Is  kept  on  a 

reduced  load  and  presumably  has  a failure  rate  smaller  than  an  opera- 
tional or  hot  standby  unit.  A cold  standby  unit  Is  assumed  not  to  lose 
operational  ability  while  on  standby.  Furthermore,  standbys  are  con- 
sidered as  being  repairable  or  not  repairable.  A repairable  unit  which 
has  failed  can  be  repaired  in  accordance  with  some  repair  time  distri- 
bution and  placed  back  on  standby  whereas  a non-repalrable  unit,  once 

It  has  failed.  Is  lost  to  all  future  system  function. 

Besides  these  various  standby  configurations  there  are  the  two 
basic  component  structures,  series  and  parallel,  whose  reliability 
provide  lower  and  upper  bounds  on  R(t)  for  a fixed  number  of  components 
(Blmbaum,  [1961]).  In  a series  structure  all  components  must  operate 
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for  successful  system  operation,  whereas  for  a parallel  structure  all 
but  one  component  can  fall  eind  still  have  system  operation. 

Caver  [1964]  considers  two  dissimilar  units  operating  redundantly 
with  constant  hazard  rate  and  arbitrary  repair  time  distributions  and 
finds  the  mean  time  to  system  failure  (MTSF) . A more  general  situation 
is  solved  by  Gnedenko  [1966]  where  n units  are  operating  In  series  with 
only  one  standby  unit.  Again,  a constant  hazard  rate  Is  assumed  for 
the  operating  unit  and  a general  repair  time  distribution  Is  assumed. 

In  this  work,  as  well  as  Caver's,  the  Laplace  transform  of  the  time  to 
failure  distribution  Is  found  but  not  Inverted.  Instead,  the  transform 
is  used  to  find  the  MTSF. 

In  Gnedenko  et  al.  [1969]  a cold  standby  system  with  n dissimilar 
components  Is  discussed  and  a sequential  method  for  finding  R(t)  Is 
shown.  This  is  further  generalized  to  allow  for  failure  properly  to 
switch  to  the  standby  unit.  Varma  [1972]  uses  similar  methods  and 
considers  n components  In  series  with  2 components  In  repairable 
standby.  Various  repair  disciplines  are  considered  with  the  hazard 
rates  constant  and  the  repair  time  of  standby  units  assumed  to  be 
arbitrary.  Kodama  [1974]  considers  a system  with  two  dissimilar  units 
having  Erlong-f allure  distributions  and  arbitrary  repair  distributions. 
Nakagawa  and  Osakl  [1974,  1975]  study  the  stochastic  behavior  of  two 
dissimilar  units  assuming  both  failure  and  repair  are  arbitrary 
functions  of  time. 
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4.2  The  n-Compartment  Catenary  System 

Consider  the  n -compartment  catenary  system  Illustrated  In 

Figure  12.  This  can  be  seen  to  represent  a cold  standby  redundant 

system  In  which  the  function  of  system  operation  Is  located  in  one  and 

only  one  compartment  at  a given  time.  Assume  the  system  is  functioning 

in  the  compartment  labeled  n at  t ■■  0 but  has  a hazard  rate  X , (t) 

n— 1 ,n 

representing  the  intensity  of  failure  and  successful  switching  to  a 
standby  component.  The  compartment  to  which  the  operating  function 
transfers  is  labeled  n-1.  Assume  also  that  this  component  can  fall  in 

such  a way  (perhaps  a switching  failure)  that  a standby  component 

» 

cannot  assume  the  system  function,  hence  there  is  system  failure.  In 
this  case  the  operating  function  has  transferred  to  compartment  zero. 
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The  hazard  rate  for  this  type  of  failure  is  denoted  1-  (t).  Similar 

un 

failure  interpretations  are  assumed  for  the  other  components. 

It  is  important  to  remember  that  the  this  compartmental 

approach  have  a fundamentally  different  meaning  from  those  hazard  rates 
normally  found  in  the  literature.  These  differences  are  seen  when  one 
compares  equations  (28)  and  (29) . Normally  hazard  rate  is  considered 
to  be  an  intrinsic  property  of  a unit  which  is  a function  of  the  time  of 
unit  operation.  The  usual  example  of  this  is  failure  due  to  wearout. 
Here  the  origin  is  the  point  in  time  at  which  the  unit  is  activated.  On 
the  other  hand,  in  the  compartmental  context,  hazard  rate  is  an 
extrinsic  property  of  a unit  determined  by  the  environment  in  which  the 
unit  is  operating.  For  example  this  may  be  a heat,  altitude  or 
radiation  gradient.  Now  the  origin  may  be  the  point  in  time  at  which 
the  system  rather  than  the  standby  unit  is  activated. 

In  reality,  a unit's  hazard  rate  is  Influenced  by  both  its 
intrinsic  properties  and  the  environmental  conditions.  If  the  environ- 
mental effect  is  assumed  to  be  constant,  this  component  of  the  hazard 
rate  has  the  "lack  of  memory"  property  and  is  therefore  insensitive  to 
the  origin  of  the  time  scale.  In  the  same  way,  if  the  unit's  intrinsic 
effect  on  the  hazard  rate  is  assumed  to  be  constant,  this  contribution 
to  the  unit's  reliability  is  Insensitive  to  the  location  of  the  origin. 
The  latter  is  assumed  to  be  the  case  in  applying  this  compartmental 
approach. 


I 

I 


.,1 


I 


I 
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4.2.1  The  Cumulant  Generating  Fxmctlon 

Using  the  cumulant  generating  function  for  the  mixed  model  (20) 
and  realizing  the  6^(t)  are  identically  zero,  one  can  write  the 
cumulant  generating  function  for  this  catenary  system 


where 


and 


and 


K(0,t) 


” 0 2 0 
ln[l+e  ^-l)p  . (t)],  tn[l+  Z (e  ^-l)p.,(t) ] , . , 
^ 1-1 


n 0, 

ln[l+  Z (e  -Dp  (t)] 
1-1  ^ 


(30) 


Pll(t)  - exp{-g^(t,0) }i  1 - l,...,n, 
'l-l  *^1+1 


1-1 


1-1 


1 * 2g3gs«.tn  and  i ■ 1,2».«.,1-1, 


X 

gj^(x.y)  - / XQj^(z)dz, 

y 

X 

g^(x,y)  - / [Xqj^(z)+Xj^_j^  j^(z)  )dz 

1 — 2,3,... ,n. 


IJ 


(Zj^,  . . . , 


^1,1+l^V 


f 
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However,  In  the  present  context  this  cumulant  generating  function 
can  be  greatly  simplified.  Recall  that  an  Initial  distribution  of 
malts  X^(0)  ■ X^,  for  all  1,  Implies 


km  - X,0,+...+X  0 . 

11  n n 


Hence  the  functional  form  of  k,  obtained  from  (31),  can  be  used  In 

equation  (30).  Also,  since  all  of  the  system  function  Is  located  In 

compartment  n at  t ■■  0,  the  Initial  distribution  Is  X (t)  ■■  X and 

n n 

X^(0)  - 0,  lihi. 

Therefore,  the  cumulant  generating  function  becomes 


n 0 

K(0,t)  - X^ln[l+^E^(e  M)Pj^(t)] 


where 


t t t 

Pln^t)  “ / / •••/  r (t  ,...t  )exp{-g  (t,t  ) 

^ n-1  '^1+1 

n-1  n-1 

- ^ 8k(tv  ® 1-1, 2,..., n-1, 

k-1+1  k-1*^ 


p^(t)  - exp{-g^(t,0)>. 


4.2.2  The  Stochastic  Distribution  of  (X, (t),  X„(t),...,X  (t)) 

1 z n 


The  cumulant  generating  function  In  (32)  Is  readily  seen  to  be 
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associated  with  an  n-varlate  multinomial  distribution  with  parameters 

p.  (t) , 1^1 •••••n  and  X . 
in  n 

Thus  the  complete  distribution  of  the  stochastic  vector 
(X^(t),  X^Ct),  ...,X^(t))  Is  determined  and  can  be  expressed; 


Prob  {(Xj^(t),  X2(t),...,X^(t))}  - 


(33) 


where  the  defined  in  (32).  It  follows  that  the  marginals 

are  also  multinomials  whose  parameters  can  be  seen  to  be  sums  of 
appropriate  probabilities,  P^(*^)  ^q* 

In  particular,  the  marginal  distribution  of  X^(t)  is  binomial  with 
parameters  X^  and  Pj^(^)*  means,  variances  and  covariances  can  be 
written  directly  from  the  parameters  In  (33)  or  as  special  cases  of 
the  expressions  In  (21) . These  are 

E(X^(t»  - X^p^(t)  , 

n\W)  - 

and 

cov(Xj^(t),Xj(t))  - -*nPin(c)Pj^(t) . (34) 
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Another  random  variable  of  particular  Interest  Is 


X._(t)  - I X (t). 
1-1  ^ 


If  one  assumes  this  system  Is  turned  on  and  operated  until  failure 
occurs  and  that  this  Is  repeated  times,  then  X^(t)  Is  the  nvmber  of 
repetitions  In  which  the  system  was  still  operating  at  time  t. 

This  variable  Is  clearly  distributed  as  a binomial  with  parameters 


This  fact  Is  exploited  In  the  next  subsection. 


4.2.3  The  Probability  of  Faultless  Operation;  R(t) 


It  was  assumed  that  X^  units  were  placed  Into  operation  at  time 
t - 0.  The  expected  proportion  of  these  still  operating  at  time  t is 
the  reliability.  Hence,  from  (34),  (35)  and  (36) 


E(X_(t))  - E(  E X.(t)) 
1-1 


n 

I EX  ft) 
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Therefore 


■ X E 
n in 


R(t)  - E 


34.(0' 


- ^ Pln<*^>  * 
1-1  “ 


(37) 


where  P^jj(0  Is  defined  In  (32). 

An  alternative  formulation  for  this  application  Is  to  consider 

to  be  an  Indicator  variable  having  the  value  0 or  1.  At  time  t - 0, 

X - 1,  then  X.  ■ 1 provided  the  system  Is  functioning  in  compartment  1 
n 1 

and  X^  - 0 otherwise.  Now  the  expected  value  of  X^(t) , and  hence  R(t), 
will  still  be  given  by  (37) . 


4.2.4  Lower  Confidence  Limit  for  the  System  Reliability  and  the 
Variance  of  R(t) 


X»(0 

An  estimate  of  R(t)  Is  given  by  R(t)  - — = 


n 

Since  the  distribution  of  X.^(t)  Is  known  to  be  binomial  It  Is  a 
standard  procedure  (see  Ostle  [1972])  to  write  the  lower  100a  per  cent 
confidence  limit  for  R(t)  as  the  smallest  value  of  p satisfying 


X \ . X -1 

[p]  [1-p]  “ ^ 1-0. 


(38) 
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Letting  this  value  be  p*  one  has 


Prob  {R(t)  i p'}  ^ 1-a. 


The  variance  of  R(t)  can  be  derived  by  considering  the  variance 
of  X,j,(t).  That  Is 


Hence 


/N 

Var(R(t))  - Var 


-^ar(X^(t)) 


4.2.5  The  Mean  Time  to  System  Failure  (MTSF) 


♦ e 

Assuming  f(t)  Is  the  failure  density  for  the  system,  then  by  the 
definition  of  expectation. 


MTSF  - E(t)  - / tf(t)dt. 
0 
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so  substituting  this  expression  Into  (40)  yields 


KTSF  - -/  tdR(t). 
0 


Integrating  by  parts  one  obtains 


KTSF  - / R(t)dt. 
0 


Hence  combining  (37)  and  (41)  one  has 


(41) 


n * 

KTSF  - I / P.„(t)dt  (42) 

1-1  0 “ 


where  , l“l*2,...,n  are  defined  In  equation  (32). 

4.2.6  The  Variance  of  the  Llfelength  of  the  System 
Using  the  definition  of  variance  one  h<w 

Var(t)  - E(t-E(t))^ 

- E(t^)  - [E(t)]^ 

00 

- / t^f(t)dt  - tE(t)]^ 

0 
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2/  tR(t)dt  - [E(t)]‘ 
0 


(43) 


Hence  from  (37),  (42)  and  (43), 


n * 

Var(t)  ••  2 E / tp  (t)dt  + 
1-10 


(44) 


where  p.  (t)  1-1, 2,..., n are  defined  In  (32). 

In 


; 


4.2.7  The  Probability  Fewer  than  r units  have  Failed;  R(r,t) 

Of  the  units  placed  Into  operation  at  time  t - 0,  the  expected 
proportion  of  these  still  operating  In  compartment  n,n-l, . . . ,n-r+l  Is 
the  reliability  of  the  first  r units  used.  This  will  be  denoted  R(r,t) 
Therefore  consider 


E(X^(t)+X^  ,(t)+...+X„  _,(t)) 
n n-i  n— r+l 


E E X (t) 
1-n-r+l 


X E P.„(t) 
" 1-n-r+l 


Hence 


E X (t) 
1-n-r+l 


R(r,t)  - E 
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i-n-r+1  (45) 

where  the  are  the  same  as  In  (32). 

Confidence  limits  on  this  can  be  derived  in  a manner  exactly  the 
same  as  the  method  used  in  Subsection  4.2.3. 

4.2.8  A Numerical  Exan^>le  Using  the  Catenary  Model 

A numerical  example  will  now  be  considered  using  a three-compart- 
ment system.  The  Joint  cumulant  generating  function  is  obtained  from 
(32).  This  is 


K(^,t)  - X )ln[l+  E (e  -l)p,,(t)], 
l-l 


where 


(46) 


t t t t 

Pl3<t)  - / / A^2^V^23^*'2^®’‘P^“/ 

2 ^2 

t- 

-/  [Xq2(z)+X22(z)  ]dz}dtj^dt2, 

t t 

P23(t)  - / X23(t2)exp{-/  [Xq2(z)+Xj^2^^^ 

° *^2 


[Xq^(z)+X22(z) ldz}dt2 


I 
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-.1 


I 

« 

I 


and 


PssCt) 


t 


exp{-/  [AQ^(z)+A23(z)]dz}, 


Now  let 

a(t)  = 

b(t)  » AQ2(t)  “ 

and 

a(t)+b(t)  » XQj^(t) . 


This  represents  a situation  in  which  each  unit  is  vulnerable  to 
two  classes  of  hazards.  Those  hazards  which  cause  failures  for  which  a 
standby  unit  may  be  provided  have  hazard  function  a(t) . Those  hazards 
without  standby  redundancy,  such  as  perhaps  hazards  to  sensing  or 
switching,  have  hazard  function  b(t).  The  last  unit  to  operate  is 
vulnerable  to  both  Icinds  of  failure  without  standby  redundancy. 

Now  making  these  substitutions  into  (46) , and  assuming 


and 


t 

A(t)  ■ / a(z)dz 
0 


t 

B(t)  / b(z)dz, 
0 


1 
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I 

i 

I 

I 

I 


I 

I 

I 


I 


one  obtains. 


p,_(t)  = / J a(t  )a(t,)exp{-/  [a(z)+b(z) ]dz}dt,dt„, 
0 0 ^ 


which  can  be  written. 


Pl3(t)  = exp{-A(t)-B(t) 


l[A(t)]^ 

. 


Similarly,  one  can  show 


p23(t)  = exp{-A(t)-B(t) }A(t) 


and 


P33<t)  = exp{-Au)-B(t)  }. 


Now  the  reliability  of  the  system  can  be  obtained  from  (37) . 
Therefore 


3 

R(t)  « £ 

1=1 


= exp{-A(t)-B(t)  } 


l+A(t)-4^|i^ 


fl 


(47) 


Now  suppose,  for  further  illustration,  a(t)  is  approximately  a 
normal  density  function  (of  course  a hazard  function  need  not  be  a 
probability  density  function)  with  a mean  of  5 ar-'  a variance  of  1. 

f a(0)  Is  essentially  zero.  The  other  hazard  function  b(t)  will  be 
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I 

! 


taken  to  be  a constant  .05.  Therefore 


and 


a(t) 


b(t)  = .05. 


One  can  now  write 

2 

Pl^Ct)  = exp{-$(t-5)-.05t}  1 

P2^(t)  = exp{-*(t-5)-.05t}  [$(t-5)]  , 
PssCt)  » exp{-*(t-5)-.05t}  » 


(48) 


where  4>  Is  the  cumulative  distribution  function  of  a normal  (5,1) . 

Figure  13  shows  a graph  of  the  probabilities  in  (48)  and  of  R(t) 
given  by  (47) . It  is  evident  that  most  of  the  failure  in  the  system 
occurs  from  t = 3 to  t « 7 when  the  a(t)  hazard  function  is  dominant. 
After  t o 7 the  constant  hazard  function  b(t)  is  causing  an  exponential 
failure  probability  to  occur. 

The  MTSF  is  obtained  from  equation  (42^  where 

3 * 

m:SF  - E / P.,(t)dt 

1-1  0 

- 18.8. 


5 


Units  of  time,  t 
Figure  13 

The  Probabilities,  PjL3(t) , the  System  Reliability  R(t)  and  the 
li.’i/.ard  FmictlouH  a(t)  and  b(t)  for  the  Catenary  Model 


73 


The  variance  of  the  lifelength  is  given  by  equation  (44)  . One 
obtains 


Hence 


Var(t)  " 2/  tR(t)dt-[MTSF]^ 
0 

2(368. 79)-(18. 80)^ 
« 384.14. 


i'Var(t) 


19.60. 


Using  the  result  in  Subsection  4.2.7  one  can  express  the  probabil- 
ity fewer  than  r units  have  failed,  R(r,t).  For  example,  the  probabil- 
ity fewer  than  2 units  have  failed  is 


3 

R(2,t)  =.  Z p (t). 
1“2  ^ 


The  expression  R(r,t)  is  simply  the  reliability  of  the  first  r 
units  of  the  system  and  many  properties  can  be  derived  concerning  this 
subsystem  by  reapplying  the  results  of  Section  4.2.  The  graphs  of 
R(l,t),  R(2,t)  and  R(t)  are  shown  in  Figure  14.  Also  indicated  on  this 
figure  are  the  MTSF  for  the  first  component,  the  first  and  second 
component  and  the  entire  system.  The  results  depicted  in  Figures  13 
and  14  can  be  found  in  Table  1. 


Table  1 


Three-Coopartment  Example  of  a Standby  System  Using  the  Catenary  Model 


t 

P33(t) 

P23(t) 

Pl3(t) 

R(l,t) 

R(2,t) 

R(t) 

0.0 

1.00 

0.00 

0.00 

1.00 

1.00 

1.00 

1.0 

0.95 

0.00 

0.00 

0.95 

0.95 

0.95 

2.0 

0.90 

0.00 

0.00 

0.90 

0.90 

0.90 

3.0 

0.84 

0.02 

0.00 

0.84 

0.86 

0.86 

4.0 

0.70 

0.11 

0.01 

0.70 

0.81 

0.82 

5.0 

0.47 

0.24 

0.06 

0.47 

0.71 

0.77 

6.0 

0.32 

0.27 

0.11 

0.32 

0.59 

0.70 

7.0 

0.27 

0.26 

0.13 

0.27 

0.53 

0.66 

8.0 

0.25 

0.25 

0.12 

0.25 

0.50 

0.62 

9.0 

0.24 

0.23 

0.12 

0.24 

0.47 

0.59 

10.0 

0.22 

0.22 

0.11 

0.22 

0.44 

0.55 

4.3  The  Mixed  Model 

Turning  now  to  the  general  mixed  model  depicted  In  Figure  7,  this 
can  be  Interpreted  as  a Competlng-Rlsk  Model  with  standby  redundancy. 
Mann  et  al.  [1974]  give  a brief  review  of  the  formulation  of  the 
Competlng-Rlsk  Model  and  reference  researchers  who  have  Investigated  It. 
It  Is  assumed  that  a unit  has  n modes  or  risks  of  failure.  One  of 
these  has  no  standby  unit  and  hence  may  represent  a complete  failure  of 
the  malfunction  sensing  or  switching  system.  The  remaining  n-1  risks 
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have  standby  units  available  to  provide  continuous  operation  of  the 
system.  It  Is  assumed  the  hazard  rate  associated  with  the  1^^  risk 
at  time  t Is  given  by  . This  time-dependent  hazard  rate  could 

describe  an  Intrinsic  risk  to  the  Initial  operating  unit  such  as 
wearout,  an  extrinsic  environmental  risk  or  some  combination  of  the 
two.  However,  for  subsequent  units  the  risk  from  intrinsic  causes 
must  be  assumed  constant. 

The  general  solution  to  this  model  is  given  In  Section  2.2.  Note 

that  the  standby  units  can  be  assumed  to  have  a further  standby 

structure  using  the  catenary  portion  of  the  model.  The  immigration 

portion  of  the  model  will  be  assumed  to  be  zero.  However,  systems 

themselves  could  be  connected  In  a probabilistic  manner  using  transfers 

such  as  X _(t). 
nO 

The  three-compartment  system  depicted  In  Figure  15  will  now  be 
developed  in  the  reliability  context.  These  results  can  easily  be 
extended  to  the  n-compartment  model.  In  the  three-compartment  case  the 
Initial  operating  unit  can  fall  In  two  ways  such  that  standby  units  are 
available  and  one  way  such  that  no  standby  Is  available  and  hence 
system  failure  occurs.  The  usual  example  In  which  two  failures  may 
occur  Is  an  electronic  device  which  falls  due  to  an  open  circuit  or  due 
to  a short  circuit.  In  general  one  could  have  type  A,  type  B or 
switching  failure.  It  is  assumed  type  B failure  can  occur  to  either  the 
initial  unit  in  compartment  three  or  to  the  standby  unit  operating  in 
compartment  two  following  a type  A failure. 


I 


I 
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Figure  IS 

Three-Compartment  Competlng-Rlek  Model  with  Standby  Redundancy 
4.3.1  The  Cumulant  Generating  Function 

Just  as  In  Section  4.2.1,  the  cumulant  generating  function  can  be 
found  from  the  n-compartment  result  In  equation  (20)  or  by  using  the 
three-compartment  result  given  by  equation  (12) . The  Initial  distribu- 
tion Is  Xj(0)  ■ X^,  Xj^(O)  ■ X2(0)  ■ 0,  and  the  ^j,(t)  are  Identically 


zero. 


*1  I 


Hence, 


^0^,02,63)  •=  X303 


gives  the  form  of  the  function  k. 

The  cumulant  generating  function  is  therefore 
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3 0. 


K(0,,02,e-,t)  - X in[l+  I (e  ^-l)p  (t)] 

^ O ij 


where 


(49) 


t r t t -N 

Pl3<t)  - J X3^3(t2>exp|-/  (XQ3(z)+Xj^3(z)+X23(z))dz-/  XQj^(z)dzUt. 


rr  C ^2 

(XQ3(z)+Xj^3(z)+X23(z))dz 
■^Q  ^ 2 3 ^ 2^  ^ ^ ^ ^0  3 ^ ■’■^3  ^ 2 3 ^ ^ * 


P23(t) 


-/  (Xq2(z)+^j^2^*^^^V*^^2  ’ 


P33(t) 


0 


exp^-/  (XQ3(z)+X^3(z)+X23(z))d^  . 
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4.3.2  The  Stochastic  Distribution  of  (Xj^(t) , X2(t) , X^Ct)) 


One  would  expect  this  vector  to  have  a quadranomlal  distribution 
and  Indeed  the  cumulant  generating  function  (49)  confirms  this  and 
provides  the  parameters  X^,  Pj^^(t),  P22(l)  and  Pj^Ct)  . Thus 

Prob  {(Xj^(t),  X2(t),  X3(t))}  - 


3 


where  the  are  defined  in  (49) . 

The  marginal  distributions  are  completely  analogous  to  the 
marginals  associated  with  the  catenary  model  discussed  In  Section  4.2.2. 
The  means,  variances  and  covariances  are  written  as  In  equation  (34). 
However  the  parameters  are  given  in  (49) . 

4.3.3  The  Probability  of  Faultless  Operation;  A(t) 

Recalling  that  X^(0)  ■ and  X3,(t)  is  defined  to  be 


X3.(t)  - Xj^(t)  + X2(t)  + X3(t), 
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the  probability  of  faultless  operation  up  to  time  t is  given  by 


R(t)  = E 


Pl3(t)  + P23(t)  + P33(t) 


where  the  are  defined  in  (49). 


(51) 


4.3.4  Miscellaneous  Properties  of  Interest 

Knowing  the  complete  multivariate  distribution  of 
(Xj^(t),  X2(t) , X3(t))  given  in  (50),  one  can  derive  many  interesting 
properties  of  the  reliability  system.  Some  of  these  were  discussed  in 
detail  in  Section  4.2  and  these  derivations  will  not  be  repeated  now. 
However,  the  final  expressions  are  identical  in  most  cases  with  only 
the  parameters,  p^j(t)  differing.  Thus  a discussion  of  confidence 
limits  can  be  found  in  Section  4.2.4.  The  mean  time  to  system  failure 
and  its  variance  are  found  in  Section  4.2.5  and  Section  4.2.6.  Other 
reliability  analyses  of  interest  may  depend  on  the  particular  problem. 

4.3.5  A Numerical  Example  Using  the  Mixed  Model 

The  three-compartment  mixed  model  illustrated  in  Figure  15  will  be 
used  in  this  example.  The  cumulant  generating  function  with  the  as- 
sociated parameters  p^3(t),  1*1, 2, 3 are  given  by  (49).  If  the  hazard 
functions  are  assumed  to  be  Identical  for  type  A,  type  B and  switching 


failures  then  these  can  be  rewritten 


r 

1 
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j i 

a(t)  = ’ 

1 

b(t)  - ” Xj^3(t)  » 

1 

c(t)  = ^03^^^  ’ 

1 

a(t)  + c(t)  » ^02^^^  ’ 

1 

1 

a(t)  + b(t)  + c(t)  ■=  XQj^(t). 

1 

1 

Suppose  also  that 

, I 

t 

A(t)  « / a(z)dz  , 
0 

1 

1 

1 

t 

B(t)  » J b(z)dz 
0 

1 

1 

and 

t 

C(t)  » / c(z)dz. 
0 

1 

Now  the  parameters  In  the  cumulant  generating  function  can  be 

1 

revnritten 

1 

P33(t)  - exp{-A(t)-B(t)-C(t)}  » 

■ 

P23(t)  ■ exp{-A(t)-B(t)-C(t) }A(t) 

1 

1 

and 

p--(t)  - exp{-A(t)-B(t)-C(t) }|A(t)B(t)-/  a(z)B(z)dzl . 

L 0 J (52) 

1 

1 

I,/  1 

f 

Assuming  now  that  a(t)  and  b(t)  are  hazards  of  a periodic  nature, 
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these  can  be  expressed  in  general  as 


and 


a(t)  = a^  + a^  sin(a^t  + a^) 


b(t)  = bj^  + b^  sln(b^t  + b^) , 


It  Is  easy  to  show 


and 


A(t)  = a^^t  + -^(cos  a^  - cosCa^t  + a^)) 


^2 

B(t)  • b^t  + ^(cos  b^  - cos(b^t  + b^)) 


Now  keeping  c(t)  completely  general  one  can  write 


p^2^t)  = exp{-C(t)  - (a^^  + b^)t  - -^(cos  a^  - cosCa^t  + a^)) 


- -r^Ccas  b,  - cos(b„t  + 
b^  4 J 

a, 

Pz3(t)  =•  p^j(t)[aj^t  + —(cos  a^  - cosCa^t  + a^))]  , 


Pl3(t)  - P2j(t)I(aj^t  + — (cos  a^  - cos(a^t  + a^)))x 

(bit  + :g^(cos  b^  -cos(bit  + b^)))-/  (Si  + a2Sin(a^z 


(b.z  + ■r^(co8  b -cos(b,z  + b, )))dz], 
i 4 J 4 


■>•  a^))x 


(53) 
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I 

1 


L 1 ^ 


For  purposes  of  this  example  these  expressions  will  be  greatly 
simplified  by  assuming  a = b,  ■=  1/5,  i = 1,2,  a =■  b.  =>  1,  a,  = 0 and 

11  3 j 4 

b^  = tt/2.  In  addltlor  let  c(t)  = 0. 

The  hazard  functions  are  therefore 

a(t)  " b , 

b(t)  « i + y cos  t I 
c(t)  = 0. 

With  these  simplified  functions,  equations  (53)  become 

PjsCt)  » exp{(cos  t - sin  t - 2t  - l)/5>  » 
p23(t)  = p^2^b)(t  + 1 - cos  t)/5  » 

r.2 

Pl^(t)  * p^^Ct)  "2  “ ^ t + t sin  t 

(54) 

These  functions  are  graphed  in  Figure  16  along  with  the  type  A 
and  type  B hazard  functions.  One  can  note  the  correspondences  between 
the  hazard  rates  and  the  probability  of  being  in  a particular  compart- 
ment. 

The  reliability  of  the  system,  R(t) , can  now  be  calculated  using 
(51) . This  function  is  graphed  in  Figure  17  and  clearly  displays  the 
cyclical  structure  assumed  for  this  model.  In  addition  using 


sin  t cos 


25 


AIR  FORCE  INST  OF  TECH  WRI6HT-PATTERS0N  AFB  OHIO  F/0  12/1 

ON  THE  distribution  OF  SOME  STOCHASTIC  COMPARTMEnTAL  MODELS  HAV— ETC(U) 
AU6  76  J 0 EPPERSON 


Units  of  tlas,  t 


Figure  16 

The  Probabilities,  Pj^3(t)  the  Hassrd  Functions  s(t)  end  b(t)  for 


the  Mixed  Model 
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00 


MTSP  - / R(t)dt 
0 


one  obtain*  a MTSF  of  approxlMtely  3.9  unit*  of  tljie.  Eaploylng  the 
reeult  given  In  (43)  for  the  variance  of  the  llfelength  of  the  system 
one  obtains 

flO 

Var(t)  - 2/  t R(t)dt  -[E(t)]^ 

0 

- 2 X 17.48  - (3.89)^ 

- 19.82 


The  standard  deviation  of  the  llfelength  Is  therefore 
*Var(t)  ■ 4.45. 

Table  2 contains  the  calculated  values  for  P23(t) > P23(^) » P33(t) 


and  R(t). 


r'*-- 


I 


Table  2 

Three-Conpartaant  Ezaa|>le  of  a Standby  Systea  Ualng  the  Mixed  Model 


t 

P33<t) 

P23<t) 

Pl3(t) 

R(t) 

0.0 

1.0000 

0.0000 

0.0000 

1.0000 

1.0 

0.5167 

0.1509 

0.0239 

0.6914 

2.0 

0.2822 

0.1928 

0.0405 

0.5156 

3.0 

0.1967 

0.1963 

0.0354 

0.4284 

4.0 

0.1687 

0.1908 

0.0342 

0.3938 

5.0 

0.1421 

0.1624 

0.0547 

0.3592 

6.0 

0.0952 

0.1150 

0.0739 

0.2841 

7.0 

0.0508 

0.0736 

0.0652 

0.1895 

8.0 

0.0266 

0.0486 

0.0456 

0.1208 

9.0 

0.0172 

0.0375 

0.0323 

0.0869 

10.0 

0.0141 

0.0335 

0.0269 

0.0745 

U.O 

0.0123 

0.0295 

0.0265 

0.0683 

12.0 

0.0089 

0.0216 

0.0254 

0.0559 

13.0 

0.0050 

0.0130 

0.0192 

0.0372 

14.0 

0.0026 

0.0076 

0.0120 

0.0222 

15.0 

0.0015 

0.0051 

0.0079 

0.0145 

• I,  * 

I :r 


I 


“ -<S!:V 
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4.4  The  General  Irreversible  Model 

The  use  of  stochastic  conpar mental  Modeling  In  reliability  can  be 
extended  to  a large  class  of  redtaidant  systea  structures  using  the 
general  Irreversible  Model  derived  In  Section  3.  Any  structure  of 
units  or  systeas  of  units  In  which  the  fmctlonal  flow  Is  Irreversible 
can  be  Modeled  In  this  way  provided  the  assun^tlons  of  Section  4.1  are 
satisfied. 

One  can  see  In  Figure  9 there  Is  a hierarchical  failure  structure 
In  the  general  aodel.  This  structure  can  be  used  to  evaluate  a Mission 
reliability  In  which  certain  coaponent  failures  or  Malfunctions  can 
only  occur  in  a specific  order.  For  exaaple,  In  analysing  the  reliabil- 
ity of  a space  craft  one  could  consider  an  Initial  class  of  Malfunctions 
associated  with  the  launch  phase.  Malfunctions  of  this  kind  nay  be 
overcone  leading  eventually  to  Mission  success;  the  Malfunction  could 
lead  to  a problea  of  another  class  during  a later  phase  of  the  Mission 
or  Mission  failure  could  occur.  In  a slallar  Manner,  If  the  Mission  Is 
undergoing  a aalfmctlon  In  a later  phase,  failure  could  occur  as  a 
consequence;  the  Mission  could  sujnrlve  to  yet  another  class  of  Mal- 
functions or  the  Mission  could  reach  untlaate  success. 

4.4.1  A Hmarlcal  Kvaaple  Pslnk  the  General  Model 

In  order  to  Illustrate  a Mission  reliability  application,  a situa- 
tion will  be  considered  wherein  diere  are  diree  distinct  phases  of  the 
Mission  with  regard  to  the  type  of  hasard  encountered.  These  phases  are 


I 


represented  by  conpartaents  1,  2 and  3 in  Figure  18. 


Trouble- 

Free 

Mission 

coiBpartaent  4 


- .05t 


^2^(0  - .05t 


.(t)  ■ .It 


A, ,(t)  - .02t 


Phase  I 
Malfunctions 

X23(t)-.06t 

P ^ 

Phase  II 
Malfunctions 

Xj^2(O".04t 

conpartiMnt  3 

T 

conpartnant  2 

Phase  III 


conparfent  1 


Xg^Ct)  - .02t 


Xn,(t)  - .02t 


Xjj^Ct)  - .004t 


Figure  18 

CoapardMntal  Eapresantatlon  of  Mission  Kallablllty  Bxa^>le 


CospartMot  4 repraascits  a troubla-fraa  slaslon.  The  only  way  In  which 
■Isslon  failure  can  occur  la  throu^  a phase  I,  II  or  III  self unction. 
The  hasard  functions  In  this  esaaple  wore  selected  prlaarlly  for 
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ease  of  computation.  However  they  were  also  selected  such  that 
phase  I malfunctions  are  more  likely  early  In  the  mission  and  phase  II 
and  111  malfunctions  become  more  likely  as  t increases.  Also  If  a 
phase  1 malfunction  has  occurred,  then  a phase  11  or  phase  111  mal- 
function Is  more  likely  to  occur. 

As  in  the  previous  examples,  the  reliability  of  this  mission  at 
time  t Is  given  by  the  probability  of  being  In  the  compartmental 
system  at  time  t.  This  can  be  determined  using  the  parameters  of  the 
ciamilant  generating  function  given  In  Section  3.2.  Equation  (25)  Is 
the  cumulant  generating  function  for  the  four-compartment  general 
Irreversible  model.  Hence  using  p^(t),  P2^(t)  and  Pj^^(t) 

from  (25)  and  the  particular  hazard  rates  given  In  Figure  18,  one 
obtains 


p„(t) 


P34(t)  • 


,-.ost^^-.uS 


PjjCt) 


2 2 
1.4286  + 0.0714  e“*^*^ 


- 1.5000  , 


Pj^(t)  - 17.9789  1.0194 


— 05t^  — 

+ 1.4881  a - 18.4476  e • 


Values  for  tltcse  probabilities  have  been  tabulated  In  Table  3. 


1 


• » 


I 


91 


Table  3 

Hlasloa  Bellabtllty  Bxaaple  Ualng  the  Geoersl  Hodel 
t **34^^^  **14^*^ 


1 

0.0 

1.0000 

0.0000 

0.0000 

0.0000 

1.0000 

1.0 

0.9048 

0.0464 

0.0241 

0.0242 

0.9996 

2.0 

0.6703 

0.1484 

0.0868 

0.0890 

0.9945 

1 

3.0 

0.4066 

0.2311 

0.1631 

0.1762 

0.9769 

1 

4.0 

0.2019 

0.2474 

0.2244 

0.2690 

0.9428 

5.0 

0.0821 

0.2044 

0.2509 

0.3574 

0.8948 

s 

6.0 

0.0273 

0.1380 

0.2391 

0.4356 

0.8400 

1 

7.0 

0.0074 

0.0788 

0.1996 

0.4995 

0.7853 

1 

8.0 

0.0017 

0.0391 

0.1484 

0.5458 

0.7350 

1 

9.0 

0.0003 

0.0171 

0.0997 

0.5734 

0.6905 

1 

10.0 

0.0000 

0.0067 

0.0610 

0.5835 

0.6513 

• 

1 

i 

11.0 

0.0000 

0.0024 

0.0343 

0.5795 

0.6162 

i 1 

12.0 

0.0000 

0.0007 

0.0179 

0.5651 

0.5837 

t 

1 

13.0 

0.0000 

0.0002 

0.0087 

0.5438 

0.5527 

1 

14.0 

0.0000 

0.0001 

0.0039 

0.5184 

0.5224 

j;  1 

15.0 

0.0000 

0.0000 

0.0017 

0.4907 

0.4923 

; ' '■ 

20.0 

0.0000 

0.0000 

0.0000 

0.3466 

0.3466 

1 

25.0 

0.0000 

0.0000 

0.0000 

0.2210 

0.2210 

1 

30.0 

0.0000 

0.0000 

0.0000 

0.1275 

0.1275 

V i 

35.0 

0.0000 

0.0000 

0.0000 

0.0666 

0.0666 

1 

40.0 

0.0000 

0.0000 

0.0000 

0.0314 

0.0314 
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In  addition,  Che  reliability  was  calculated  using 

A 

R(t)  - 2 

i-1 

and  diese  values  are  Included  on  Che  table. 

Mow  one  can  find  the  aean  tiae  to  aission  failure  by  applying 
equation  (41).  This  results  in 


/ R(t)dt  “ 10.8  mlts  of  tiae. 
0 


The  variance  of  the  expected  aission  lengtii  is  provided  by  equation 
(43) . Hence 


and 


Var(t)  - 2/  tR(t)  - (E(t)l^ 
0 

- 282.56 


»Var(t)  ■ 16.8. 

Slallarly,  these  principles  can  be  applied  to  the  probability  of 
being  in  the  failure-free  state  at  tiae  t,  to  get  the  aean  and 

variance  of  the  tiae  to  first  aalfunction.  Thus 
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and 


E(tlae  to  first  aalfunction)  - / p .(t)dt 

0 

“ 2.8  units  of  time 


Var(tlme  to  first  malfunction)  ■ 2/  t p..(t)dt  - [2.8]^ 

0 ^ 

- 2.16. 


The  values  from  Table  3 have  been  graphed  In  Figure  19.  One  can 
see  the  relatively  short  time  of  expected  operation  without  malfimctlon 
by  the  small  area  tnder  the  P^^(t)  curve. 

Also  from  time  0 to  4,  If  the  system  has  malfunctioned.  It  Is 
slightly  more  likely  to  be  a phase  I malfunction,  but  from  time  4 
onward  the  phase  III  malfunctions  becoM  much  more  likely.  Beyond 
about  15  time  units,  either  the  mission  has  failed  or  the  system  Is  In 
a phase  III  malfunction  state. 


Probability,  a Mission  Is  Malfunction-free  or  In  Malfunction  Phase  I,  II  or  III 
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5,  CONCLUSIONS 


5 , 1 Sumnary 

The  cumulant  generating  functions  are  derived  for  a compartmental 
model  which  Is  a mixture  of  the  catenary  and  mansnlllary  models  axid  a 
general  Irreversible  model.  The  first  and  second  moments  are  found 
under  general  assumptions  concerning  the  initial  cumulant  generating 
functions  and  In  particular  when  pulse  labeling  Is  assumed  the  exact 
multivariate  distributions  are  Identified. 

Substructures  of  these  models  are  applied  to  several  reliability 
situations.  A restricted  form  of  the  catenary  model  Is  used  to  analyze 
a standby  redundant  system.  The  mixed  model  Is  used  for  a redundant 
system  in  which  there  are  multiple  risks  of  component  failure.  And  the 
general  model  Is  used  for  a mission  reliability  analysis. 

Several  fundamental  properties  of  Interest  in  reliability  are 
derived  In  terms  of  the  hazard  rates.  Some  of  these  are  the  system 
reliability,  confidence  limits  of  the  reliability,  the  mean  time  to 
system  failure  and  the  variance  of  the  llfelength  of  a system.  A 
numerical  example  Is  given  of  each  of  these  configurations  to  Illustrate 
some  of  the  concepts  developed. 
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5.2  Future  Research 

5.2.1  Other  Considerations  Conceralag  Rellafelllty 

In  Section  A,  a number  of  questions  arise  which  are  worthy  of 
future  study.  By  what  experimental  design  could  one  best  determine 
the  hazard  functions  of  a given  system?  It  Is  likely  cost  would  be 
an  Important  factor  In  collecting  these  data,  hence  optimal  designs 
In  terms  of  time,  or  equipment  would  be  of  Interest. 

Can  one  develop  a compartmental  approach  In  which  the  hazard 
function  contains  both  an  age-'dependent  component  and  a time-dependent 
component? 

Cardenas  and  Matls  [1974]  derive  a two-compartment  reversible 
model  with  time-dependent  transition  rates.  A study  should  be  made  to 
consider  the  use  of  this  model  In  reliability  where  unit  repair  Is 
considered. 

5.2.2  Considerations  Concemlng  Organizational  Manpower 

The  problem  of  labor  turnover  has  received  attention  from  several 
researchers  and  many  methods  for  labor  turnover  analyses  have  been 
proposed.  Lane  and  Andrew  [1955]  develop  a survival  curve  method  where- 
in a group  of  new  employees  are  observed  until  departure  from  the 
organization.  Inferences  are  then  made  based  on  the  type  of  curve 
obtained  In  this  way.  Herbst  [1963]  Is  mentioned  In  Section  1 with 
regard  to  his  successful  use  of  a deterministic  compartmental  approach 


I 

I 

I 
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In  studying  the  labor  turnover  phenomenon.  Clowes  [1972]  simplifies 
the  Herbst  model  by  using  fewer  compartments.  His  fit  of  actual  data 
is  still  good  and  his  approach  Is  also  deterministic.  Several  other 
models  for  social  systems  can  be  fotmd  in  Bartholomew  [1973].  More 
insight  into  the  labor  turnover  process  might  be  gained  by  using  a 
stochastic  compartmental  approach.  It  would  be  particularly  interesting 
to  use  this  approach  to  study  the  labor  turnover  process  in  a graded 
social  structinre  - such  as  the  military  rank  system. 


5.2.3  Considerations  ConceminK  Survival  Distributions 


The  subject  of  survival  distributions  Is  closely  related  to  the 
field  of  reliability.  Indeed,  one  of  the  most  recent  works  on  the 
subject  Is  by  Gross  and  Clark  [1975]  who  take  the  approach  that  theory 
developed  to  study  the  life  of  a mechanism  or  a system  can  In  many 
cases  apply  Just  as  well  to  the  life  of  an  animal  or  a human  being. 

There  are  two  concepts  developed  in  Section  4 which  could  be  pursued 
In  the  area  of  survival  distributions. 

Firstly,  the  catenary  model  could  be  used  to  model  a multlstaged 
disease  in  which  death  occurs  with  differing  probability  in  each  stage 
or  transfer  to  a more  advanced  stage  can  occur.  Such  a model  would  be 
particularly  well  suited  in  instances  where  an  outside  influence  affects 
the  transition  rate  from  stage  to  stage. 

An  example  of  this  Is  a seasonal  effect  or  some  other  periodic 
Influence  on  the  advancement  of  a disease.  The  existence  of  biological 
rhythms  Is  well  established  (see  Metis  at  al.  [1973])  and  the 
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advancement  of  some  Illnesses  could  be  related  to  these  variables 
using  a compartmental  model. 

Secondly,  the  mamillary  model  could  be  used  to  Investigate 
competing  risks.  The  theory  of  competing  risks  Is  well  advanced 
(see  Chlang  [1968])  and  It  may  well  be  that  a stochastic  compartmental 
approach  would  produce  little  or  no  new  theory.  However,  the  fact 
that  the  compartmental  method  Is  different  may  allow  one  to  obtain 
previous  results  in  a simpler  manner  or  provide  Illumination  one  might 


not  gain  otherwise. 
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APPENDIX  A:  DERIVATION  OP  THE  COMDLAHT  GENERATING  FUNCTION 

FOR  THE  GENERAL  lOISL 


Th«  aolutton  to  (22)  wtU  now  bn  dnrlvwl.  Thn  nubaldlnry  equations 
can  be  written 


1 -e,  " T3  

(l-e  )XQ^(t)  ^ b^jj(t)+(l-e  ^)Xjjj(t) 

dx  

“ Q g f J “ 2,3, ..«,n  . 

E (e  ^-l)A.o<t) 

J-1  (Al) 


Now  letting  - (e  ^-1),  1 - 1,2,..., n.  It  follows  that  dV^  - e ^d0^. 
Using  this  relationship  and  the  subsidiary  equations  (Al),  the  following 
systea  of  ordinary  differential  equations  In  the  and  K can  be 
obtained.  This  systea  Is 


•*^1 " '^iN)i<‘>‘*^ » 

"1  - ♦ *t\,i<‘>| 


dtt  i-2,3,...,n. 


dX(e,t)  - ZV.X,n(t)dt. 
J-1  J 


(A2) 


Now  solving  (A2)  sequentially  for  the  and  K(6,t)  In  tenv  of  the 
trensltloo  rates  ^^j(t)  end  nfl  constants  of  Integration,  this  systea 
can  be  solved  for  the  constants  of  Integration  to  obtain 
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Cj^  - Vj^  eatp{-hj^(t,0))  » 
t 

C2  ■ **Pt“\C*:,tj^)-h2(tj,0)}dtj^4V2  expi-hjCt.O) } » 


‘3  ■ 

J-1  t t 
W 0 

2t  t t 

/ / •••/  (t^  J^) 

°Vl  S-1  ^ ^ 

•-2 

1-1  *i*^i+i  *1+1  * Vi*^  J * * *1  ^ 


m-2 


■>1 


“ J*  h»  — I’^i.  » fO) J n dt|  1 ^ •*• 

1-1  1 1 * *1+1  * J 1-1  *1  ^ 

t t t 

+/  / •••/  ^IJ  ^*1**2*  * * * * *J“1^ 

° *j-l  *1+1 


J-1 

- £ t 
t-2 


J-1 


] 


J-2 
+ E V 
1-2 


"N^*i-i**j-i^“'*j^*j-i»°^^^*i-i^*j-i  ■*■  ••* 
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and 


+ y t t t 


»-2 


1-1  ^I’Vl  Vl  ^ Vl’^  J-1  **'  1 * 1^-1" 

■“2  B_1 

1-1  1 \ ^ 1+1  ^ J J ^ 1-1  ^ 


t t t 

^1J^^1’*^1+1 Vl^  exp{-h^(t,t^) 

” J-1  Vl 


j-1  J-1 


-1  - 

1-1+1  * J J-1  1 


J " • • • fO» 


+ Vj  e*p{-hj(t,0)}  f 


''nfl 


- K(l.t)  - Vjj/%^Q(t^)  exp{-h^U,t^)>dt^ 

\ilVt  ^iO<V^ll^‘l-l> 


n A 

— <-S<*"*k-l>-V*k-l>*t>-*«+l<'.-'n>  - 


+ 

I<Ai<A2< 


O t t t 

2 III 

.T?<A ,<n  0 t t-  , 

" Vr‘ 


"^'W:r 
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m-2 


J'^1  1 1-1  *1+1  * * 

m-2  B-1 

1-1  *1  *1  * *1+1  * Vi  Vi  ^ " 1-1  “ 

1 1 t 

^Vt  ’ “^t  ^n0^‘tt^*^ln^*l»*2 Vl^ 

n 2 

“p<-'‘i<'-'i>-°“i+i<‘i-'i+i>>j«i]-  - 

- V2[^„V2.0<V  “P<-V2<‘>V>'*'» 

*11  Vl,0<W2.«-l<V2‘  “f‘-V2<'"'n-2> 

^*1 


t t 


-Vi<'..-2-‘.»-'V2‘'vn  »„o(v».^2.«<w 

n 

«xpt-h  .(t,t  ,)-h  (t  ))<It  ,dt 

n-z  n-1  n n-1  n n-1  n 


c t t 


^Vt  ^nO^‘“^V2,n-l^*n-2^Vl.n^*ii-l^ 

o n— 1 


“P«-V2<'-'n-2>-Vl<V2-'.-l>-‘‘n<*.>-l'V> 


dt  .dt  .dt 
n-2  Q-1  n 


rt 


•] 


■ Vl<*>  /p^n-l,0<*n>  •*P^-Vl<**VJ‘**« 


t t 


*IJ^  ‘oo‘W2.a^Vl> 

II 


"^<*«t-l**a»‘‘*n 


'-■1' : 
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where 


and 


n 0 no  n a n n 


X 1-1 

h (x,y)  - / E X,-(x)d* 
y J-0 


^ij^*l*“***j-l^  ■ \,l+l^V**’^j-l,j^*j-l^‘ 

6l 

Now  substitute  back  e -1  for  and  let  «nd 

P^^(t)  - ezp{-h^(t,0) } I 1 - l,2,...,n, 

t 

Pl.1+1^^^  “ -^O^.l+l^V  ®*P^-*»i(‘.V"*‘l+l^*'l»0>^‘*’=l»  ^‘1»2 n-1, 

‘'w+lVtj  exp{-h^(t.t^^) 

M t t t 

'*’  / / / •••/  (tg  j) 

Kt  <i,<.V.<l  <J  0 t,  , t.  , ^*1 


1 2 •••  «-l  •*  ''  "j-1  -l^-l 


•”2 

1-1  *1**1+1  *1+1  * Vl’J  J ^ * *1  ^ 

*-2  ^1 
1-1  *1  *1  * *1+1  * J J ^ 1-1  ^1  ^ 


t t t 

"'L  *’ij^‘i*‘i+i Vi^  **p<-hi(t.ti) 

j-i  *1+1 


ll; 


I 
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fi  (t)  - / X «(t  )exp{-h  (t,t  )}dt  , 
n Qnun  n n n 


5„  ,(t)  - / A (t  ) exp{-h  -(t,t  )}dt  +/  / ^_n<t  )X  . .) 

Q^l  Q H“ljO  tl  TX^X  Q ^ 0 t ® ^ XfQ  H“X 

n 


*n-1.0*V^n-2,n-l*'ii-2*“’’*"''n-2''’'«-2' 

n 


Vz’  V V ^n-Z,n^ Vl^ 

t t t 

+/  / / „ ,(t  JX  - „(t  J 

Of  ^ nO  n n-2,n-l  n-2  n-l,n  n-i 

*^n  'n-1 

Vz  ' Vz^  Vz » Vl^  "*‘n^  Vl  ’ V Vz*^  Vl**^ 


and  tn  genaral 


n t t 


•»I’<-*j<‘"'l-l*"^t('l-l’'n)’'‘tl-l’l'n 
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n-J  J+Jl-1  t t t 

***’{“hj  (•=•  t^)  > 

■‘’  * * * 

+ S /'/'  /'  •••/'  h “<V*1  « -i’ 


m-2 


ii*2 


” ^9  9 _i)«*p{-h  (t,t,  _,)-  E h (t,  .,t  .; 

1-1  *i’*l+l  \+l  J *-1  J-  i.i  \ 1 1 

m-1 

"^1  ,“l’^n\"/*^)l  -l**^n  + ... 

»-l  B-1  i-1  1 

'"K  K ^nO^V^Jn^*^j’^j+l**”' Vl^ 

“ n ^j+2  *^J+1 

n-1  n 


) 


exp{-hj  (t,  tj)-^E^hi^j^(ti,  ti^j^)  }^n^dti. 


(A4) 


Mow  (aJ)  can  be  written  more  concisely  as 


Cl  - Uj^Ci.t)  - (e  -l)pj^j^(t)  » 

6 9 


j ®1 

Cj  " “jCfi.t)  • E c«  -l)p<<(t). 


1-1 


U' 


n e. 


Ci^l  - “ K(e,t)-^E^(e  ^-l)«j(t) 


(AS) 


> f 


I 
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.r-’- 

I 


where  the  P^j(t)  and  ^j(t)  are  given  by  (AA) . 

by  characterlatlc  theory,  these  n^l  constants  can  be  related  by  a 
functional  relationship  i|«  such  that 


(A6) 


the  form  of  iji  will  be  deteroilned  by  the  Initial  conditions.  When  t - 0 
the  functions  u^  becoae 

64 

''‘4(1*0)  « e -1  l«l,2,...,n 

and 

u„+lCi.O)  - K(1*0)  - k(^  . 

Hence  when  t ••  0,  equation  (A6)  reduces  to 

9 9 9 

kC9)  - *(e  ^-1,  e ^-1 e “-!)  . (A7) 

O4 

In  order  to  write  <l>  In  the  usual  functional  form  let  y^  « a -1.  This 
Implies  that  9^  <■  In(y^-fl)  and  (A7)  now  becomes 


>^(^“(74+1) *^n(yj^+l))  - ♦(y4,y2....7j^).  (A8) 


The  function  i>  has  now  been  defined  by  the  Initial  form  of  the 
Joint  cumulant  generating  function.  Therefore  from  equations  (A6)  and 
[A8] 


Ill 


Uj^ld.t)  - k(ta(uj^d.t)+l),  ....  £n(u^(0,t)+l)), 


The  solution  to  (22)  can  now  be  obtained  by  substituting  (AS)  Into  (A9) 
to  obtain  (23) . 


I 
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APPENDIX  B:  DERIVATION  OP  THE  HDMENTS  FOR  THE  MIXED  MODEL 

AND  GENERAL  MODEL 

The  BOMnts  will  be  derived  by  differentiating  (20)  or  (23)  with 
respect  to  the  appropriate  paraaeters  and  evaluating  this  for  9^-0 
i " 1, . . . ,n.  Define 

^ ®i 

Uj  - in[l+^2:^(e  -l)Pj^j(t)]  for  J - 1,2,..., n. 

Therefore  (20)  can  be  written 

(Bl) 

for  1 - l,2,...,n. 

(B2) 

The  second  partial  derivative  of  (Bl)  with  respect  to  is 


(B3) 


Hence 


n 0. 

K(i,t)  - h(u)  + r (e  J-l)«,(t) 

J-1  ^ 


I 
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I 


Similarly  one  obtains 


36^36^^ 


All  e'B  - 0 


-Pik(t)Pj^(t)  IJ  5 k,  l^J 

0 1 or  J > k , 


Aasumlng  the  means,  variances  and  covariances  for  the  Initial 
distribution  are  as  defined  In  Subsections  2.1.3,  2.2.3  and  3.1.4,  one 
can  cosine  (B2),  (B3)  and  (B4)  with  (B5),  (B6)  and  (B7)  to  obtain  the 
expressions  for  E(X^(t)),  V(X^(t)  and  cov(X^(t) , Xj(t))  as  shown  In  (21) 
and  (24). 


115 


VITA 


Jon  Ohmn  Epperaon  waa  bom  In  Alva,  Oklahom,  Septeaber  16,  1937. 

Ha  la  the  aon  of  Raverend  Mltchall  Stokaa  Bpparaon  and  Mra.  Eathar 
Rookar  Epparaon.  Kra.  Epperaon,  now  widowed,  llvea  In  Corralaa,  Mew 
Hexlco.  Her  nailing  addreaa,  and  permanent  addraaa  for  Jon  Epparaon, 
la  P.O.  Box  574,  Corralaa,  New  Mexico,  87048. 

Before  Jon  waa  one  year  old,  Mr.  Epperaon  waa  called  to  the 
Flrat  Freabyterlan  Church  in  Ada,  Oklahona  where  Jon  attended  public 
achool  until  the  age  of  12.  In  January,  1949,  the  Epperaona  moved  to 
Albuquerque,  New  Mexico  and  Jon  graduated  fron  Albuquerque  High  School 
In  1956.  For  two  yeara  Jon  attended  The  College  of  Wooater  at  Wooater, 
Ohio,  then  tranaf erred  to  die  Unlveralty  of  New  Mexico  at  Albuquerque 
where  he  conpleted  a Bachelor  of  Science  degree,  with  a najor  In 
natheaatlca,  in  January,  1962. 

Jon  entered  the  U.S.  Air  Force  in  March,  1962  and  haa  alnce  been 
atatloned  at  Little  Rock,  Arkanaaa;  Coluahia,  Mlaaourl;  Haablngton,  D.C.; 
Colorado  Springe,  Colorado  and  College  Station,  Texaa.  At  the  Unlveralty 
of  Mlaaourl,  he  conpleted  a Maater  of  Arta  degree  in  nathenatlca  in 
January,  1968.  After  three  and  a half  yeara  at  the  Pentagon  and  two  on 
the  nathenatlca  faculty  at  the  U.S.  Air  Force  Acadeaqr,  Jon  cane  to 
Texaa  AAM  Unlveralty  to  aeek  a Doctor  of  Fhlloaophy  degree  in  the  Inatl- 
tute  of  Statlatlca.  He  expecta  to  complete  thla  degree  by  Auguat,  1976 
and  then  return  to  the  nathenatlca  departnant  at  the  Air  Force  Acadeny. 


The  typlat  for  thla  dlaaertatlon  waa  Mra.  Suaan  Epperaon 


